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Abstract 

The behavior of solutions of the finite-genus Whitham equations for the weak dispersion Hmit of the 
defocusing nonlinear Schrodinger equation is investigated analytically and numerically for piecewise- 
constant initial data. In particular, the dynamics of constant-amplitude initial conditions with one or 
more frequency jumps (i.e., piecewise linear phase) are considered. It is shown analytically and numer- 
ically that, for finite times, regions of arbitrarily high genus can be produced; asymptotically with time, 
however, the solution can be divided into expanding regions which are either of genus-zero, genus-one 
or genus-two type, their precise arrangement depending on the specifics of the initial datum given. This 
behavior should be compared to that of the Korteweg-de Vries equation, where the solution is devided 
into the regions which are either genus-zero or genus-one asymptotically. Finally, the potential appU- 
cation of these results to the generation of short optical pulses is discussed: the method proposed takes 
advantage of nonlinear compression via appropriate frequency modulation, and allows control of both 
the pulse amplitude and its width, as well as the distance along the fiber at which the pulse is produced. 

The weak dispersion limit of the defocusing nonlinear Schrodinger (NLS) equation has been extensively 
studied in recent years (see e.g., Refs. L13t l22l l27l l29l [31 1), and in some sense it is well-characterized 
mathematically. Not as much is known, however, about the detailed behavior of the solutions for specific 
choices of initial datum I^E^tSHHl purpose of this work is to present analytical and numerical results 
regarding the behavior of a special class of solutions of the NLS equation in the weak dispersion limit. More 
precisely, we consider the initial value problem for the NLS equation in the weak dispersion limit with the 
initial data having a constant-amplitude and piecewise constant frequency. We believe that, on one hand, 
this behavior is interesting mathematically, and that, on the other hand, it could have potential applications 
in the generation of intense, ultra-short optical pulses. As such, it is worthy of further study. 

The structure of this document is as follows: in section^we introduce the problem, and in section|2lwe 
review some well-known results regarding the weak dispersion limit of the NLS equation. In section |3l we 
discuss the behavior of solutions corresponding to "single-jump" initial conditions, which are the starting 
point for our investigation. Then, in sections |4] and |5] we present the analytical calculations which are the 
main results of this work. In section |6l we demonstrate these results through numerical simulations of the 
NLS equation and obtain further information about the solution behavior, and in section we discuss the 
application of our results to the generation of intense, ultra-short optical pulses. Appendix lA. II describes 
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our nondimensionalizations and our choice of units, Appendices IA.2I and IA.3l review some known results 
regarding genus-one (i.e., periodic) solutions of the NLS equation, the Whitham averaging method and 
the NLS -Whitham equations, and Appendix IA.4I gives the details of some calculations whose results are 
presented in sections |3l and 0] 



1 The NLS equation with small dispersion 

In this section we recall some basic results regarding the behavior of solutions of the nonlinear Schrodinger 
(NLS) equation with small dispersion. This will establish the background and the notation necessary to 
extend these results in the following sections. 

The semiclassical limit of the NLS equation. We start from the defocusing NLS equation (that is, the 
NLS in the normal dispersion regime in optical fibers) in dimensionless form: 

where we assume < £ ^ 1. In the context of optical fibers, t represents the dimensionless propagation 
distance and x is the dimensionless retarded time (cf. Appendix lA.il) . To study the weak dispersion limit of 
the NLS Eq. (Il.lt . we first express the field q{x,t) in a WKB form as 



q{x,t) = a/pCm) exp[/(p(x,0/e] , (1-2) 

where p{x,t) = \q{x,t)\^ and cp(x,f) = £ arg q{x,t) = {ie/2)\n[q* {x,t)/q{x,t)] represent respectively the local 
intensity and the normalized local phase of q{x,t). We then introduce the normalized phase gradient 

u{x,t) = d(^{x,t)/dx, (1.3) 

that is, u{x,t) = {ie/2) {ql{x,t)q{x,t) — q* {x,t)q^{x,t)) /\q{x,t)\^ . Throughout this work we refer to x and 
t as the space and time variables, respectively. By analogy with the fiber optics context, however, we will 
refer to u{x,t) as the local frequency of the solution. With this decomposition, the NLS Eq. (Il.lt can be 
written in the form of a conservation law: 



dp d{pu) 



dt dx 



(1.4a) 



d(pu) dfTlnlnd^ 
^^^--lpu^ + -p^--e^p — lnp . (I.4b) 



dt dxy"^ ' 2^ 4 ^3x2 

The weak dispersion limit of the defocusing NLS equation is defined as the problem of studying the solutions 
of Eq. dl.ll) as £ ^ 0+ with initial condition expressed in terms of Eq. (11.21) . The limit is singular, and it 
should be considered in the weak sense. 

There are two main situations where the weak dispersion limit of the NLS equation is relevant: nonlin- 
ear fiber optics and Bose-Einstein Condensation (BEC). In the context of nonlinear fiber optics, the weak 
dispersion Umit is relevant for the long-distance transmission of non-return-to-zero (NRZ) pulses, as dis- 
cussed in Ref. l25l . or for the generation of intense short optical pulses, as proposed in this work. In the 
context of BEC, the semiclassical limit applies due to the very small value of Planck's constant H relative 
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to quantities associated with macroscopic objects, i.e., h = e. Normalizations appropriate for long distance 
optical fiber communications were discussed in Ref. 1 14|, whereas in appendix lA.il we discuss scalings and 
nondimensionalizations relevant for the generation of intense short optical pulses. We emphasize however 
that the results presented in this work apply equally well to Bose-Einstein condensates and to other physical 
context where the NLS equation is relevant, such as for example ferromagnetics and water waves. 

Hydrodynamic analogy and dam-breaking problem. If both p and u are smooth and p > 0, and if 8 ^ 1, 
Eqs. (II. 4t are approximated to leading order by the following reduced hydrodynamical systeml^HIHl 

dt (u) " (l m) 3^ (m) ■ ^^-^^ 

Equation (II. 5t is called the dispersionless NLS equaton, and is used to describe a surface wave motion in 
shallow water. In the hydrodynamical setting, p and —u represent respectively the depth and velocity of 
water, and x and t are dimensionless space and time. For p > 0, the eigenvalues u ± of the coefficient 
matrix are real, and the system (II. 5t is strictly hyperbolic. This system, which is known as the shallow water 
wave equations and has been intensively studied (see e.g. Ref. [341), can be rewritten in Riemann invariant 
(i.e., diagonal) form 

3- =5^:3-, k = \,2, (1.6) 
at ax 

where the Riemann invariants ri 2(x,f) are given by 

ri=u-2^/p, r2 = u + 2^, (1.7) 
and the characteristic speeds si^2{x,t) are 

■^1 = ^(3ri +r2) =u-^, S2 = ^(n +3r2) =u + ^. 

Note that Sk> implies a left-moving wave (i.e., dx/dt = —Sk), and that Eqs. M.ll are equivalent to 

P = ^(r2-ri)2, „=l(ri+r2). (1.8) 



Since the system of PDEs described by Eq. (11.61) is strictly hyperbolic, it is possible to show that, for 
"rarefaction" initial data, namely, when ri 2(^,0) are both monotonically decreasing functions of x, a global 
solution exists for all f > 0. In many cases of interest, however, the initial data do not satisfy this property 
and as a consequence they develop a shock, as we shall see in the following. 

Square-wave initial conditions. Consider first the initial datum given by the following rectangular pulse 
of width 2L: 

P(..0) = {^» (1.9a) 

m(x,0)=0, (1.9b) 

with qo > 0, corresponding in the optics framework to an NRZ pulse. The system of equations dl.St with 
initial conditions (11.91) is known in the literature as the "dam-breaking" problem. (In the hydrodynamic 
analogy, the problem describes the behavior of a mass of water which is initially confined in a uniform. 
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spatially localized state by two dams located at x = ±L and both of which are removed at ? = 0.) The 
Riemann invariants for this situation are shown as the dashed lines in Fig. 12. Ih . Note that the initial data 
for the Riemann invariants ri , r2 corresponding to the initial condition (ll.9t is not of rarefaction type. It 
is possible however to obtain initial data of rarefaction-type by properly redefining the initial value of the 
invariants r\ and ri for \x\ > L, as shown in Fig. 12. lb later. (This procedure is a special case of the process 
known as regularization, see next section.) Then the system has the following solution up to the time 
fo = L/qQ. for < ;c < L + lq^t, it is 



while for ;c > L + Iq^t, p(x,f) = u{x,t) = 0, with p(— x,f) = p(x,f) and u{—x,t) = —u{x,t) (cf. Ref. ll25l ). 
The full solution of the NLS equation (11.11 ) with initial condition (ll.9t . as obtained from numerical simu- 
lations (described in section |6ll is depicted in Fig. ll.lh . (This kind of solution is usually called ufan in the 
context of hydrodynamics.) The speeds of the boundaries of the top (p(x,f) = q^) and bottom (p(x,f) = 0) 
regions are easily obtained from Eqs. (ll.lOt (see also Fig. 12.1b later); these two speeds are respectively 
^2 = and S2 = —2qQ. 

Equations (ll.lOt cease to be valid beyond the time to = L/qo when the boundaries of the top region 
meet at x = (i.e., the time at which the left-moving characteristic emanating from x = L meets with the 
right-moving characteristic from x = — L). An analytical expression for q{x,t) when t > L/qo however can 
be obtained using the dispersionless limit of the scattering transform for the NLS equation' '. The corre- 
sponding behavior of the full solution of the NLS equation is shown in Fig. II. lb . The appearance of small 
oscillations in the numerical solution in Figs. ll.lh .b was discussed in Ref. fT4l . and is the consequence of ap- 
proximating the discontinuous initial data (ll.9t with a continuous initial datum in the numerical simulations 
(as described in section |6ll. It should also be noted that some care must be taken regarding the regulariza- 
tion of the discontinuous initial datum (ll.9t for Eq. (II. 5t . since a weak solution with a discontinuity is not 
unique. In fact, one can construct a different solution of (ll.4t with a discontinuity. In our case, however, 
the regularization described in section [SJenforces the continuity of p(x,f), thus removing the ambiguity and 
producing a unique solution. 

Frequency jumps and high-frequency oscillations. In terms of optical pulses, the above results imply 
that the initial condition in Eq. ( I1.9t rapidly spreads out, as shown in Figs. Il.lh .b. This behavior can be 
partly prevented (or, alternatively, reinforced) by employing initial conditions with nontrivial phase. For 
example, consider the following: 



with p(x,0) still given by Eq. (I1.9at . Hereafter, we will use r°(x) and r2(x) to refer to the value of the 
Riemann invariants ilJl at f = 0. If p(x,0) is given by Eq. ( I1.9at and u{x,0) by Eq. ( I1.9bll i. for |x| < L we 
have 




(I.IOa) 
(1.10b) 




x<0 
x>0 
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x<0 
x>0 



x<0 
x>0 



(1.11b) 



(1.11a) 
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as shown by dashed Unes in Figs. I3.2h .b. Note that the value of the invariants can be redefined for \x\ > L, 
since p{x,0) = there. 

In terms of optical pulses, Eq. ( ll.9bU amounts to imposing a frequency jump at the center of the pulse 
{x = 0). If Mo > 0, the right half of the pulse acquires a positive frequency and the left half a negative 
frequency. Thus, owing to the normal dispersion, the two halves of the pulse will tend to move towards each 
other. In terms of the hydrodynamical problem, this corresponds to assigning an inward initial velocity 
to the mass of water, as if two pistons were acting on each side of it. (For this reason, this case is often 
referred to as the "piston" problem.) Note that if mq > the initial data r*' 2(jc) are increasing. Thus, another 
consequence of the initial frequency jump is that if mq > a shock develops at ;c = 0, and the solution 
develops high-frequency oscillations, as shown in Fig. 13. Ih . (This type of shock is called collisionless, or 
dispersive, to distinguish it from the usual type of shock, which is dissipative; e.g., see Ref. |34|.) The 
characteristic frequency of these oscillations is 0{\/z) (i.e., one period of the oscillation shown in Fig. 13. lb 
is of order s). If mq < instead, the two halves of the pulse will move away from each other, and no shocks 
develop in this case. (In hydrodynamics, solutions such as this one are called of rarefaction type.) When 
the solution develops a shock, the hyperbolic system (II. 5t ceases to be valid, and the solution of the NLS 
equation in the weak dispersion limit must be obtained by properly regularizing the hyperbolic system, as 
we briefly discuss next. 

2 Regularization and the NLS-Whitham equations 

The semiclassical limit of the NLS equation has been extensively studied in the last fifteen years; e.g., see 
Refs. [I3l |22l |23l 123 Ell ISl ESl El for different approaches, the connection with the integrable character 
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Figure 1.1: Deformation of an NRZ pulse (ie., the "dam-breaking" problem) defined by the initial conditions in 
Eqs. (|r9}, with = 1, L = 20 and = 0.1: (a, left) t' = 10; (b, right) t' ^ 40, where t' = tlz. The critical time 
here is fo = L/qo = 20; hence Fig. II. Ih and Fig. II. lb show the solution respectively before and after the time fo- 
The dotted lines show the initial condition, while the solid lines show the result of numerical simulations of the NLS 
equation ( II. Il l, which are performed in terms of the fast time scale f ', as discussed in section|6l 
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and the multi-phase solutions of the NLS equation and with Whitham's averaging method. A self-contained 
description of the regularization process can be found in Ref. |25|, together with a detailed treatment of 
the situation in which the initial condition contains only one frequency jump. Here we will limit ourselves 
to presenting a brief general overview and recalling some results which are relevant for the remainder of 
this work. Some additional details (which are necessary to perform the calculations described in sections |3l 
andlU are contained in Appendices IA.2l and lA.3l 

As we will see throughout this work, for certain kinds of initial conditions the solution of the NLS equa- 
tion (ll.lt with small dispersion develops high-frequency oscillations. When this happens, the approximate 
hyperbolic system (II. 5t is inadequate to describe the dynamics of the solution, because of strong dispersive 
effects appearing due to the presence of high-frequency oscillations. An effective way to describe the be- 
havior of the solution of the NLS equation in these situations is obtained by taking an average over these 
high-frequency oscillation via the Whitham technique, which consists in locally approximating the the solu- 
tion by finite genus solutions of the NLS equation, then describing the global behavior as a slow modulation 
of this local periodic or quasi-periodic structure. The evolution of these modulations is then governed by 
the Whitham equations, which are obtained by averaging the conservation laws of the NLS equation over 
one period of the fast oscillations, and which express the local average of the quasi-periodic solutions with 
respect to these fast oscillations (cf. Refs. [ 12 J3^^| for the Korteweg-de Vries equation). 

Recall that the hyperbolic system (ll.6t is a dispersionless limit of the NLS equation. The presence of 
high-frequency oscillations in the solution of NLS corresponds to the formation of a shock singularity in 
Eq. (11.61) generating strong dispersion, and the Whitham averaging technique then provides an appropriate 
dispersive regularization of this singularity. This regularization turns out to consist in enlarging the sys- 
tem il.6l to include more than two Riemann invariants (as given by Eqs. (I2.3t below) in such a way that 
the initial data for the enlarged system becomes of rarefaction type, which in turn implies that the system 
possesses a global solution for all values of time. The solution of the original problem is then described in 
terms of the solution of an NLS-Whitham equation of finite genus. The dispersionless equation ( ll.6t corre- 
sponds to the genus-zero NLS-Whitham equation. The value of the genus is determined by the specifics of 
the initial condition considered, since these determine the number of Riemann invariants which are needed 
to regularize the hyperbolic system. 

Let us briefly review some features of the regularization and introduce the NLS-Whitham equations (see 
Appendices IA.2l and lA31 and Refs. fT3||22l|25l|22l for more details). It is well-known that the NLS equation 
is integrable via the inverse scattering transform. The scattering problem associated with the defocusing 
NLS equation (ll.lt is given by the eigenvalue problem ElUU l \ = z\, where v is a two-component vector, 
z G C is the spectral parameter of the scattering problem, and where the Lax operator L is defined by 

A genus-^ solution of the NLS equation is associated to a genus-g hyperelliptic Riemann surface R:w^ = 
figiz), with 

2g+2 

k=i 

The solution of the NLS equation dl.lt corresponding to Eq. (12. 2t is described in terms of a Baker- Akhiezer 
function constructed from Riemann theta functions with g phases (e.g., see Ref. ||5lEl for details). The 
branch points ri , . . . ,r2g+2 of R determine the spectrum of L. Since the Lax operator L is self-adjoint, these 
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branch points are all real, and we label them so that ri < r2 < • • • < ?'2g+2- The spectrum of L is then given 

by (-00, n] U [r2,r3] U • • • U [r2g,r2g+i] U [^2^+2,°°). 

For an exact genus-^ solution of the NLS equation, the branch points rj , . . . , r2g+i are obviously con- 
stant, independent of space and time, owing to the isospectrality of the inverse scattering transform of the 
NLS equation. The solution then gives a quasi-periodic solution with g phases of the NLS equation, and 
with a small dispersion of order £^ implies that each period of the phase is of order £ (i.e., high-frequency 
oscillations of order 1/s). 

Suppose now that the ^-phase solution is slowly modulated. Then the spectral parameters r]^ are expected 
to shift slightly (and even additional small gaps are expected to open in general). The shifted parameters are 
constant as the eigenvalues of the Lax operator L . We can treat the modulation problem using a singular 
perturbation method as follows: if we consider modulations on the scale of (x^t) in the NLS equation (ll.lt . 
the modulations are of order one, and the oscillations are of order 1 /s. This implies that one can treat 
the modulation problem as a small perturbation of the g-phase solution, which is the key of the Whitham 
averaging method. This method is based on an adiabatic assumption that the leading-order solution of 
the NLS equation (i.e., the ^-phase solution) is preserved under the modulation, but with slowly changing 
parameters. Employing a standard averaging perturbation method, one introduces the fast and slow time 
and space scales respectively as [y! = x/e, t' = t/e) and {x,t). Then the first-order correction to the g-phase 
solution in the perturbation method compensates the motion of the parameters so that the parameters appear 
to be constant over regions of order one in {x,t), but at the same time acquire small constant shifts due to the 
modulations in general. The corresponding equations for the spectral parameters with respect to the slow 
scales {x,t) are called the NLS-Whitham equations, or simply the Whitham equations, and are obtained by 
averaging the conservation laws of NLS with respect to the fast oscillations at order £ (see Appendix lA. 3l for 
more details). 

It was shown in Refs. |[T3ll27l that the genus-g Whitham equations for the NLS equation can be written 
in the Riemann invariants form, in which the spectral parameters give the Riemann invariants. 



for = 1 , . . . , 2^ + 2. It then follows that for each £ 7^ the solution of the NLS-Whitham equations describes 
the slow modulation of finite-genus solutions of the NLS equation. It was also shown (cf . Lemma 4. 1 in 
Ref. (25\ ) the following, which is the most important property of the NLS-Whitham equations: 

Proposition 2.1 The characteristic velocities s^ possess a double sorting property, i.e., Vk, I = 1 , . . . , 2g + 2 

they satisfy the two conditions 



Then, as a consequence of 12. II we have (Corollary 4.2 in Ref. igU): 

Corollary 2.2 If the initial values of the r/^ are each nonincreasing and if they satisfy the separability con- 
dition 




(2.3) 




(2.4b) 



(2.4a) 



(2.5) 
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\/k = \, . . . ,2g + \ at t = 0, the initial data is of rarefaction type, and therefore the hyperbolic system of 
equations (I2.3t has a global solution for t >0, i.e., it is regular 

In general, the initial conditions for the two Riemann invariants obtained from the dispersionless sys- 
tem (ll.4t do not satisfy monotonicity and the separability condition (I2.5t . In other words, Eqs. (I1.7t (that 
is, the system (I2.3t with ^ = 0) in general do not have a global solution. The regularization process then 
consists in enlarging the set of Riemann invariants so that the resulting NLS-Whitham equations have a 
global solution. This is done by representing the initial data for the NLS Eq. (Il.lt in such a way that all 
the Riemann invariants ri , . . . , ^2^+2 are monotonically decreasing functions of a; at f = 0. Note that it is 
always possible to do so for piecewise-constant initial data, since in this case there is some ambiguity in 
how the spectrum of the Lax operator is represented in terms of the Riemann invariants. This ambiguity 
can then be exploited to redefine the initial datum for the Whitham Eqs. (I2.3t by adding degenerate gaps 
(e.g., see Figs. l3.2l and l4.lt . All of the solutions discussed in sections|3j|l]and|5lfall within the framework of 
piecewise-constant initial conditions. In this case the data at ? = are always genus-0, but highly degenerate, 
in the sense that they can be described in terms of a higher genus with degenerate gaps. 

It should be noted that the adiabatic approximation (the use of the Whitham equations for the slow 
evolution) implies that the local genus of the solution is preserved. A separate issue is how the genus of the 
solution changes from one region to the next. This is a bifurcation problem through a critical point, and can 
be approached by regularization, i.e., by trying to patch two Whitham systems with different genus at this 
point. Our regularization acts as a "globalization", in the sense that the Whitham equations with a proper (in 
general larger) genus now describes a global behavior beyond the perturbation range. For fixed x, a regular 
scheme must solve a connection problem in t because of the change of genus, and for fixed t one also needs 
to solve a connection problem to match regions with different genus. 

It is also important to realize that for piecewise-constant initial data there is more than one way to 




Figure 2.1: Evolution of the Riemann invariants in two equivalent cases. Dashed lines: the Riemann invariants at 
f = 0; solid lines: the invariants at f 7^ 0. Figure lzTk (left) corresponds to the square-wave initial datum in Eqs. M.9\ 
and Fig. 11.11 regularized by genus- 1 data; Figure IzTb (right) shows an equivalent diagram of the left one, and it is 
given by genus-0 data. Hereafter, the subscripts "— " and "+" refer to the value of the invariants respectively to the 
left and to the right of their initial discontinuity. The ellipses in Fig. l2.1h indicate "locking" points, i.e., points for 
which r\ (x*,f) — r2(jc*,f) — r^,{x^,,t) (on the left) and r2(ji;*,f ) = r^{xif,t) — r4(xH.,f) (on the right) for all t. This implies 
that the regularization at those points is trivial. Finally, note that the initial conditions r^jC-^) in the equivalent figure 
(Fig. 12. lb ) have been redefined whenever p(x,0) = (i.e., u{x) itself is not defined in the NLS-Whitham equation 
when p(x) — 0, see Eqs. ( IA.9» . 
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regularize the initial datum (e.g., see Figs. I2.1h .b). The minimum number of Riemann invariants that are 
necessary so that the system becomes regular is related to the genus of the solution of the NLS-Whitham 
equations. That is, 2g + 2 invariants correspond to a genus-g solution of Eqs. (I2.3t . The local genus of 
the solution of the NLS equation is roughly speaking the number of distinct frequencies which are locally 
present in the solution over regions of order one in {x,t). Since the Riemann invariants are the branch points 
of the spectrum of the finite-genus solution of NLS which locally approximates the full solution, it is then 
clear that the local genus is equal to the number of gaps determined by the local value of the Riemann 
invariants (e.g., see Figs. 15. H and 15. 2t . The opening or closing of one of the gaps for some values of {x,t) 
corresponds to a local change of genus in the solution of the NLS equation. Note however that all finite- 
genus solutions of the NLS equation with non-zero genus become singular in the limit £ ^ 0+ . In this sense, 
the solution of the Whitham equations represents a weak limit, since when g ^0 the solutions of the NLS 
equation only converge in an average sense (i.e., weak convergence). 

Finally, with regards to Fig. I2.1l we should note that genus- 1 data is necessary in order to preserve the 
value of u{x,0). With step initial data, however, no gap opens during propagation, which means that the data 
is degenerate, producing a genus-0 solution. The situation would be different in the case of non-step initial 
data (e.g., if the transition from p(x,0) = to p(x,0) = were continuous); in that case oscillations would 
appear, as described in Ref. I.14J . 

3 Single-jump initial conditions 

We now briefly summarize some results from Ref. l25ll relative to a single-jump initial datum, since they 
provide the basis for the framework that will be used to analyze the more complicated scenarios discussed in 
the remainder of this work. We will consider the constant-amplitude wave given by the (single-jump) initial 
condition in Eqs. (I1.9at and (ll.9blt . where we take L ^ oo. The value of the original Riemann invariants 
in the genus-0 system Eq. (11.61) at f = is again given by Eqs. (II. lit , which are now valid Vx G M. Four 
different situations arise depending on the size of the frequency jump 2uq, as shown in Figs. l3.1h -d: 

(i) Mo > 2^0 > (Lemma 4.3 and Theorem 4.4 in Ref. |25 1). 

Since mq > 0, the original Riemann invariants r° 2 (x) are increasing functions of x, and therefore the 
genus-0 system (ll.6t does not have a global solution. In this case the problem is regularized by 
considering the genus- 1 NLS-Whitham equations. (That is, four invariants are necessary so that the 
resulting system has a global solution.) The Riemann invariants are related to the solution at ? = as 
follows: 



with the upper/lower signs corresponding to the upper/lower inequality for x, respectively. The qual- 
itative evolution of these Riemann invariants is shown in Fig. I3.2h . The solution develops a region 
of genus- 1 high-frequency oscillations in the central portion of the pulse, surrounded by a genus-0 
region, as illustrated in Fig. l3.lh and Fig. l3.3h . As shown in Fig. l3.lh . the genus- 1 portion of the solu- 
tion describes slow modulations of high-frequency oscillations, as would be the case in a wave packet. 
The genus-1 portion of the solution is located in the region |x| < s^t. The characteristic velocities are 



n = —Mo — 2^0 Vx, 

r2 = — Mo±2^o forx>0, 



r3 = Mo ± 2^0 for x > 
r4 = Mo + 2^0 Vx , 
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^3 = — ' '^3^ ~ ' ^i'-'^ ■^s > ■^a^ > and 

= (mo + ^o) [l + Mo^o/(wo + ^o)^] , s'^ = uq[ \ - 3a{\ -a)Ki -2a^K2]/{l -aK\) (3.1) 
(cf. Eqs. (4.16) and (4.17) in Ref. |25|), where a = Iq^/uQ < 1 and 

y(l+asin2e)(l+acos2e) / " ^ (1 +a2sin2e)(l +a2cos2e) 

Hereafter, the superscripts "— " and "+" refer to the value of the Riemann invariants respectively to the 
left and to the right of the discontinuities at f = (cf. Fig. l3.2t . Note that, owing to Eq. (12.31) . positive 
values of s correspond to left-moving Riemann invariants, and that for regularized, non-increasing 
invariants (x) > r+ {x) \/x G R. 

(ii) < Mo < 2<?o (Lemma 4.5 and Theorem 4.6 in Ref. |"25l). 

Again, since uq > 0, the genus-0 Whitham equations do not have a global solution. This case is 
regularized by the genus-2 NLS-Whitham equations. The six Riemann invariants are related to the 
solution at f = as follows: 

ri = — mq — 2^0 r4 = — mo + 2^o Vx, 

r2 = ±Mo — 2qo for x^O, = ibwo + 2^o for x^O, 

ri, = Uq — Iqo Vx , r6 = Mo + 2^o Vx . 

The qualitative evolution of these Riemann invariants is shown in Fig. 13.2b . Even though ^ = 2 is 
necessary to regularize the initial data, no genus-2 region appears. The solution develops a genus-0 flat 
region (a degenerate genus-2 region) in the central portion of the pulse, surrounded by genus- 1 high- 
frequency oscillations, as illustrated in Fig. 13.1b . The genus-0 and genus- 1 portions of the solution 
are respectively located in the regions |x| < s'^t and s'^t < |x| < s^t, with = — ^2 , s'^ = —5^, and 
where > > 0. The characteristic velocities are (Eqs. (4.28) in Ref. 1125 1 ) 

s^=qQ-^UQ, s^ = {uQ + qQ)[l+UQqo/{uQ + qQf]. (3.2) 

Furthermore, in the genus-0 region, the solution takes the value P = ^^(''2 ~ ^1)^' 

p = ^g[l + Mo/2^o]^ (3.3) 

Note that the amplitude of the genus-0 (non-oscillatory) region at the center increases with increasing 
modulation strength uq. At the same time, however, the dependence of s'^ on mq implies that the width 
of the genus-0 region decreases with increasing modulation strength, and the region ceases to exist 
when Mo > 2qQ. The maximum possible amplitude that can be obtained in this way is thus p^ax = ^Iq, 
obtained for mo = 2qQ. 

This case, < mq < 2^o. is the most interesting case for applications because, unlike the high- 
frequency oscillations in the genus- 1 portion, the high-amplitude genus-0 region can survive an ap- 
propriate filtering and produce short, high-intensity optical pulses. As mentioned above, however, the 
maximum amplitude that can be obtained with this arrangement is limited, i.e., q = {I + mo/^o)?o < 
2^0- In the next section we will see how this limit on the maximum pulse amplitude can be overcome 
by employing more than one frequency jump. 
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(iii) -Iqo < Mo < (Theorem 4.7 in Ref. ||25l). 

Since mq < 0, the Riemann invariants r\ ^ (■^) decreasing functions of x, and therefore the genus- 
NLS-Whitham equations (ll.5t with the Riemann invariants defined as in Eq. (I1.7t has a global 
solution, and no regularization is necessary. The solution develops a depression zone in the center, as 
illustrated in Fig.|3lTJ;. More precisely (Eqs. (4.31) and (4.33) in Ref. L25J ). 



p(x,0 



^g[l + Mo/(2^o)]^ 
[x+(mo + 2^o)?]V9?^ 



^0 



? *\ *\ 15*2 1 5 





Figure 3.1: Deformation of an NRZ pulse in the presence of an initial frequency jump: (a, top left) uq = 2.01; (b, top 
right) Mo = 1; (c, bottom left) mq = — 1; (d, bottom right) mq = —2.5. The dotted line shows the initial condition; the 
solid line shows the result of numerical simulations of the NLS equation, performed as discussed in section|6l In all 
four cases it is = 1, ^ = 35, and =0.1, and the solution is shown at f' = 20, where f' = t/e. 
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where 

^2=1 ("0 + 2<7o) , S2 =qo-uo. (3.4) 

(iv) MO < —2^0 (Lemma 4.8 and Theorem 4.9 in Ref. ll25l '). 

This case is regularized by the genus- 1 NLS-Whitham equations (11.51 ). The Riemann invariants at 




Figure 3.2: Qualitative diagrams illustrating the evolution of the Riemann invariants: (a, left) mq > 2qo, with a single, 
expanding genus- 1 region; (b, right) < mq < 2^o, with an expanding genus-0 region surrounded by expanding genus- 
1 regions on either side. Dashed lines: the original invariants 2(x) at f = 0; solid lines: the regularized invariants 
at f 7^ 0; dot-dashed vertical lines: boundaries between regions of different genus. Note that the connecting segments 
between in Fig. l3.2h .b. in Fig. l3.2h and in Fig. 13.2b are actually curved, and are only represented here by 
straight lines for simplicity (e.g., see Ref. |7 1). 




space X space x 



Figure 3.3: The boundaries between regions of genus-0 and genus- 1 in the (x,f) -plane corresponding the cases shown 
in Fig. 13.21 (a, left) uq — 5/2, corresponding to case (i) and Figs. 13.1^ and 13. 2h : (b, right) uq — 1, corresponding to 
case (ii) and Figs. 13. lb and 13. 2b . Dashed lines: the same boundaries after adding a constant frequency offset Mavg = 3 
to both sides of the jump, as discussed in section|4] 
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f = are defined by: 

ri=UQ±2qQ forx^O, ^3 = — mq — 2^0 Vx 

r2 = mq + 2^0 Vx , r4 = — mq i 2^0 for x ^ , 

The corresponding solution is shown in Fig. I3.1H . Even though ^ = 1 is necessary to regularize the 
initial data, the gap never opens, and this case produces a degenerate genus-0 solution. The solution 
is similar to the one in the previous case, except that q{x,t) tends to zero at x = 0. More precisely, 
p(x,f) = for |x| > s^t and p(x,f) = in the limit e ^ 0+ in the region |x| < s^t, where (Eqs. (4.36) 
in Ref. f25i) 

sl = -uo-2qo, -Mo + 2^0- (3.5) 

(The nonzero value of q{0,t) in Fig. 13. IH is due to the finiteness of £ in the numerical simulations.) 
Furthermore, the solution does not develop genus- 1 regions (the high-frequency oscillations visible in 
Fig. 13. IH disappear in the limit 8 ^ 0+ in the case of step initial data). 



The calculation of the characteristic speeds for the Riemann invariants is based on the formulation of 
the NLS-Whitham equations. We refer the reader to appendices IA.3I and IA.4I and to Refs. H |25l for 
further details. Note that the outer boundaries of the genus- 1 region are given by x = it^'outer^ where Souter 
is the same in case (i) and case (ii): Souter = ("0 + ^0) [l + mo^o/("o + ^0)^] (cf. in Eq. (I3.lt and s'^ 
in Eq. (13. 2t ). Note also that the location of the boundaries between regions of genus-0 and genus- 1 in the 
numerical simulations shown in Figs. l3.1h -d agrees very well with the analytical results just presented, even 
though the value of £ used is not very small. 

Case (ii) is the most interesting for applications because of the high-amplitude genus-0 region. Also, in 
cases (iii) and (iv) (i.e., when uq < 0) the solution is only of genus-0, but lower-amplitude. For this reason 
we will mainly focus our attention to the case of positive frequency jumps. Hereafter, we refer to case (ii) 
as a "sub-critical" frequency jump and to case (i) as a "super-critical" frequency jump. More precisely. 

Definition 3.1 We say that a single frequency jump located at x = xq is super-critical if r^i^x^) > Tj^Xq ), 
sub-critical if r^{xQ) < rjix^), where as usual r^(x^) = lim\^^o± r^(x + Ax). 

In the following sections we generalize the above solutions and discuss the behavior of solutions of the 
finite-genus NLS-Whitham equations in the presence of an arbitrary number of jumps in the initial data. 



4 Behavior of finite-genus solutions: Two frequency jumps 

More complicated situations than those described in the previous section arise when the initial conditions 
contain more than one frequency jump. The main purpose of this and the following section is to study 
the interaction among finite-genus solutions of the NLS-Whitham equations generated by those frequency 
jumps. First of all, however, let us briefly discuss the effect of adding a non-zero average frequency across 
the jump. The Galilean invariance of the NLS equation implies that it is possible to redefine the local 
frequency m(x,0) up to an arbitrary additive constant. The effect of such a constant, representing a global 
frequency translation, is just a change in the overall group velocity of the solution. More precisely, if qo{x,t) 
is a solution of Eq. dl.lt . so is 

qc{x,t) = exp[/(cx-c^f/2)/£] qo{x + ct,t), (4.1) 
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for any constant c. Corresponding to Eq. (I4.1I) there exists a Galilean symmetry of the NLS-Whitham 
equations (I2.3t for the Riemann invariants: namely, if rj{x,t) is a solution of the Whitham equations (I2.3t . 
then so is the translated solution rj{x + ct,t) +c. Note from Eq. (I4.lt that a positive frequency c/e cor- 
responds to a negative shift — c in the group velocity dx/dt. Also, as an effect of the Galilean transforma- 
tion (I4.lt the rescaled frequency of the transformed solution is shifted by c, i.e., Uc{x,t) = £3;f[arg qc:{x,t)\ = 
c + £3;i[arg <7o(x + cf,f)]. In terms of the decomposition of the (jc,f)-plane into regions of different genus, we 
thus have the following: 

Lemma 4.1 The boundaries between genus-0 and genus-1 regions for a single frequency jump upon adding 
a nonzero average frequency Mavg to the initial condition are given by x = s±t, where s± = ibsoid + Mavg . cind 
where Sou is still given by Eqs. (I3.lt . (13. 2t . (13.41) and (13. 5t in cases (i-iv), respectively. 

Thus, adding a nonzero frequency offset to both sides of a single frequency jump has the effect of tilting 
the corresponding boundaries between regions of different genus, as shown in Fig. l3.3h . If only one jump is 
present, it is obviously possible to choose the average frequency so that both of the boundaries of the outer- 
most genus-0 region move in the same direction. In the presence of several frequency jumps, however, is not 
always possible to do so, as will be discussed later. Because of the Galilean invariance, we will sometimes 
describe the initial condition for u{x, 0) in terms of the frequency jumps, defined as Cj = u{X^ , 0) — u{Xf , 0) 
for all 7 = l,...,N, where A'^ is the total number of jumps and Xi,... ,Xi\i are the jump locations. 

We now turn to initial conditions with two frequency jumps. We will consider constant-amplitude initial 
conditions with |^(x,0)| = q^Mx & M. For simplicity we will take symmetric jumps of size C\ = C2 ='■ 2uo 
located at X2 = —Xi = AX/2 > 0, so that Iuq is the size of each frequency jump, as in the single-jump case. 
That is, let 

{-2uo x<-AX/2, 
\x\<AX/2, (4.2) 

2uo x>AX/2, 

again with |^(x,0)| = ^0 > Vx G M. In the absence of either one of the two jumps, the solution would 
behave according to the theory described in Section |3l with the only exception that the average frequency 
across the jump is now nonzero. More precisely, since the average frequency at x = ±AX/2 is ztwo, the 
solutions described in the previous section would move with velocity ^uq. In other words, the central 
portion of the solutions described in Section|3]will tend to move towards each other if mq > and away from 
each other if uq < 0. As we will see, this description provides an accurate picture of the overall solution at 
sufficiently small propagation times. After this initial stage, however, significant interaction effects appear. 

In analogy with the calculations described in Section |3l let us now proceed to analyzing the solution by 
looking at the Riemann invariants. The initial values 2(jc) of the non-regularized Riemann invariants for 
the genus-0 NLS-Whitham equations (ll.6t are 

{—2qo — 2uo x<—AX/2, ( 2qQ — 2uo x<—AX/2, 

-2qo |x|<AX/2, r^{x) = <.2qo |x|<AZ/2, (4.3) 

—2qo + 2uo x>AX/2, |^2^o + 2mo x>AZ/2. 

These invariants are shown as dashed lines in Figs. I4.1h -d. Since expressing the values of all the regularized 
Riemann invariants at f = in terms of the initial datum would be rather tedious, and since all the values 
can be easily inferred from Figs. I4.1h -d. we omit the formulae for brevity. 

Based on Proposition 12. H and CoroUarv 12.21 the following four scenarios arise depending on the size of 
the frequency jumps: 
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(i) uq > 2qo > 0: 

This case is regularized by the genus-2 NLS-Whitham equations (i.e., with six invariants), as il- 
lustrated in Fig. 14. Ih . where the (two) original and (six) regularized Riemann invariants are shown 
respectively by dashed and solid lines. 

The regions of genus-0, 1 and 2 for the solution of the NLS equation are shown in the bifurcation 
diagram in Fig. l4.2h . Both jumps are super-critical individually, so two genus-1 regions open up, one 
at the location of each frequency jump. These two regions interact and they create a genus-2 region at 
the center of the solution. The genus-2 region then expands forever (cf. Propositions 14.41 and I4.7t . In 
terms of the numerical simulations described in the next section, this case corresponds to Fig. 16. Ih . 

(ii) ^0 < Mo < 2^0 : 

This case is regularized by the genus-4 NLS-Whitham equations (i.e., with ten invariants), as illus- 
trated in Fig. l4. Ib .d. where the original and regularized Riemann invariants are shown at two different 
values of t. 

The regions of genus-0, 1 and 2 for the solution of the NLS equation are shown in the bifurcation 
diagram in Fig. 14.2b . Both jumps are sub-critical individually; so two genus-0 regions open up ini- 
tially at the location of each frequency jump (each surrounded by genus-1 regions, as in case (ii) of 
section |3ll- As these two portions of the solution interact, a genus-0 region forms temporarily in the 
central portion of the solution. This region disappears after a while, however, and a genus-2 region 
forms which expands forever as in the previous case (cf. Propositions 14.41 Wj\ and lOl . In terms of 
the numerical simulations, this case corresponds to Fig. 16. lb and Figs. l6.3b .d. 

(iii) < Mo < ^o: 

This case is also regularized by the genus-4 NLS-Whitham equations. The original and regularized 
Riemann invariants are shown in Fig. 14. lb . and the regions of genus-0, 1 and 2 for the solution of the 
NLS equation are shown in the bifurcation diagram in Fig. l4.2H . 

As in the previous case, the individual jumps are sub-critical, so two genus-0 regions open up initially 
at the location of each frequency jump. Also similarly to the previous case, a genus-0 region forms 
in the central portion of the solution as the two genus-0 portions of the solution interact. Contrary 
to the previous case, however, the genus-0 portion of the solution now persists and expands forever 
(cf. Propositions 14.41 l4J\ and 1481 . Numerical simulations corresponding to this case are shown in 
Fig. 16.1b and Figs. I6.3h .b. 

(iv) Uq < 0: 

As in the single-jump case, two expanding depression regions form at the location of each frequency 
jump. Eventually, these regions merge to form a unique depression zone. As in the single-jump case, 
the details of the regularization process vary depending on the value of uq. Because the solution is 
always genus-0, however, this case is not as interesting as the previous ones in terms of applications, 
and therefore it will not be investigated further. 



Let us now discuss more quantitatively the behavior of the solutions in each of the above scenarios. 
When Mo > as in cases (i)-(iii), the initial values of the non-regularized Riemann invariants r^2i^) 
increasing functions of x, and one has: 

maxr^{x) = (oo) = lim (;c) = 2(mo - <?o) , (4.4a) 

xes. 

rmnr^{x) = ^^(-oo) = lim 4{x) = -2{uo - qo) . (4.4b) 
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As we will see in section |6l cases (ii) and (iii) offer the most useful behavior for practical applications. 
In particular, the case uq = qo, which is the separatrix between cases (ii) and (iii), is especially interesting 




Figure 4.1: Qualitative diagrams illustrating the evolution of the Riemann invariants: (a, top left) uq > 2qo, corre- 
sponding to case (i) and showing expanding genus- 1 and genus-2 regions; (b, top right) < mq < qo, corresponding 
to case (iii) and showing an expanding genus-0 central region; (c, bottom left) qo < uo < 2qo corresponding to the 
intermediate case (ii), and showing the invariants at small time t and a temporary genus-0 region in the center; (d, 
bottom right) the same intermediate case (i.e., < UQ < 2^o as in (c)), but now showing the solution at large time t 
and an expanding genus-2 region in the center Dashed lines: the original invariants r'^2i^) at f = 0; solid lines: the 
regularized invariants at f 7^ 0; dot-dashed vertical lines: boundaries between regions of different genus. As in sec- 
tion|2j the ellipses in Figs. l4.lt .d indicate "locking points", i.e., locations corresponding to the same value of x Vf 
(cf. Remark l43l 5). Finally, note that, as in Fig. 13.21 those portions of the invariants that coincide in the limit e ^ 0+ 
can be omitted when calculating the local genus. 
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since it produces a high-amplitude genus-0 region of constant width at the center of the pulse, as shown in 
Fig. 14.2b . This case is regularized by the genus-3 NLS-Whitham equations, corresponding to Fig. 14. lb with 
rj and re deleted and where now rj(oo) = r2{—°°), since uq = qo- Hereafter, we will call the case uq = qo 
the "critical" two-jump case, and we will refer to cases (ii) and (iii) as the super-critical and sub-critical 
two-jump cases, respectively. More precisely. 

Definition 4.2 We say that an arbitrary collection of positive frequency jumps Ci , . . . , located at posi- 
tions X\, . . . ,Xfj is collectively super-critical ifr^l{°°) > r2{—°°), collectively sub-critical ifr^{oo) < r2{—°°), 
and collectively critical if r^{oo) = r2{—°°), where as usual r^(iboo) = \imx^±oorf(x). 

Note the difference between the above and Definition 13.11 which distinguishes whether a single jump is 
individually sub-critical or super-critical. 

It should be noted that, even though the NLS-Whitham equations with g = 3,4 are necessary to regularize 
the data in some cases, only the genus g <2 appears in the solution with two frequency jumps (similarly to 
the scenario with a single frequency jump). Note also that, unlike the case when a single phase is present, 
calculating analytically the precise location of the boundaries between different multi-phase regions after 
the individual portions of the solutions have come into contact is a highly nontrivial task, and we have 
not attempted to do so. However, it is possible to calculate the velocity of all boundaries between genus- 
and genus- 1 regions. The main reason for the difference between the two types of boundaries is that, 
upon removal of degenerate gaps, the calculation of characteristic speeds in regions of genus- 1 is expressed 
in terms of elliptic integrals and the boundary with genus-0 turns out to be a straight line (see Lemma 
l4.3l below'). whereas the characteristic speeds in regions of genus-2 require the evaluation of hyperelliptic 
integrals (e.g., see Ref. 1 17 1). 

The explicit calculation of a few characteristic speeds is reported in appendix IA.4I Since such calcula- 
tions are rather lengthy, however, and since the methods used are described in Refs. |6l0|25l, here we omit 
the details relative to all the different cases, and we limit ourselves to summarizing the main results. 

Lemma 4.3 All boundaries between genus-0 and genus-1 regions are given by straight lines. 

Proof. The result follows from Eqs. (12. 3t and by noting that the characteristic speeds are constant since all 
the invariants are constant in the genus-0 portion. □ 

This result is obviously independent of the number of frequency jumps present, and is just a special case of 
a generic feature of hyperbolic systems: if one side of the initial condition is constant, there is a solution in 
terms of simple waves, for which the characteristics are straight lines'^. 

Proposition 4.4 The boundaries between genus-0 and genus-1 regions before the interaction are given by 
straight lines with the same velocities as in Eq. (I3.lt and Eqs. (13. 2t upon application of the appropriate 
Galilean shifts. That is, Si^± = zt^,- + mq for the regions relative to the frequency jump to the left and s,- -t = 
±5; — uofor those relative to the frequency jump to the right. 

Proof. The result follows by direct computation of the characteristic speeds in the NLS-Whitham equa- 
tions (I2.3t . but is consistent with Lemma 1411 regarding the application of a constant offset to a single 
frequency jump, showing that before interaction the behavior of the solution in each region is unaffected by 
the presence of the other jump. □ 
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Remark 4.5. In case (ii) (cf. Figs. l4.1b .d and Fig. 14.2b ). it is r5{x^,t) = r(,{x.^,,t) = rj{x^,t) for all t, where 
x^{t) = max {a; G M : r5{x,t) = 2{uq — <?o)}- That is, the three invariants are "locked" together at that point. 
Similarly, if one defines x„{t) = mm{x G R : r(,{x,t) = -2(mo -<7o)}, it is r4{x„,t) = r5{x„,t) = r^{x„,t) 



for all t. This locking phenomenon is identified by the two ellipses in Figs. 14. Ib .d (see also Fig. 12. Ih ). 
Moreover, since two of the three invariants always coincide in a neighborhood of x^ or x*^,, in each case 
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Figure 4.2: Bifurcation diagrams illustrating the regions of genus-0, genus-1 and genus-2 in the (x,r)-plane for the 
cases shown in Fig. 14. II namely: (a, top left) uq — 2001/1000 (mq > ^qo), corresponding to case (i) and Fig. l4.1h : (b, 
top right) Mo = 5/4 {2qo > mq > qo), corresponding to case (ii) and Figs. l4.1t .d: (c, bottom left) the critical case uq — I, 
which is the separattix between cases (ii) and (iii); and (d, bottom right) mq = 5/8 {qo > uq > 0), corresponding to 
case (iii) and Fig. 14. lb . The horizontal axis is position, x, and the vertical axis is time, f . In all cases, qo^l. The thick 
(blue) lines indicate the boundaries between regions of genus-0 and genus-1 (cf. Propositions 14.41 147l and l4~8i . The 
thin (black) lines indicate boundaries between genus-1 and genus-2 regions, and were computed by approximating 
them with circular arcs of given starting point and given initial and final slopes. (These conditions determine the 
circular arc uniquely.) The thin dashed lines indicate boundaries between genus-1 and genus-2 regions in those cases 
when the final slope is unknown. The dot-dashed horizontal lines indicate the values of t corresponding to the diagrams 
of the Riemann invariants in Figs. l4.1h -d. 
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we locally obtain a trivial regularization. That is, near x = x^, the thi^ee Riemann invariants r^, and 
r-i are reduced to a single invariant f{x,t) given by r{x,t) = r-i{x,t) for x < x^{t) and r{x,t) = ri{x,t) for 
X > x^{t), and a similar situation arises for r4, and near x = x^^. Those trivial regularizations are 
however necessary in order for all the Riemann invariants to be nonincreasing functions in x (unlike the case 
in Fig. l2.1h . where trivial ones can be removed). 

Conjecture 4.6 The width of the outermost genus-2 regions present in cases ( ii) and ( Hi) after interaction 
tends asymptotically to zero. 

In both case (ii) and case (iii), the speeds of the left- and right-boundaries of the genus-2 region in question 
are respectively given by s^ and Sj (cf. Figs 13.21 and 14.21) . where as before the superscripts "— " and 
refer to the values of the Riemann invariants respectively to the left and to the right of their discontinuities at 
t = 0. The fact that the width decreases monotonically and tends asymptotically to a constant follows from 
the double sorting property of the characteristic velocities. It is not possible, however, to prove whether the 
asymptotic width is zero without explicitly calculating the speeds of the boundaries between genus- 1 and 
genus-2 regions, which is beyond the scope of this work. (The appearance of this phenomenon was first 
noted by F.-R. Tian for the case of the KdV-Whitham equations with two jumps'221.) 

Proposition 4.7 After interaction, the boundary between the outermost genus-1 regions and the surround- 
ing genus-0 regions is given by x= ib (xq + i'outerO' where ^outer = in case (i) and ^outer = ■^'9 in cases (ii) 
and ( iii). In all these cases, the resulting velocity is 

Soutev = (2mo + '7o) [ 1 + 2uQqo/{2uQ + qof ] . (4.5) 

The characteristic speed can be obtained from degenerate g = 2 calculations (cf. Appendix IA.4t . Note 
however that the appropriate values of and however can also be obtained from s'^ in Eq. (I3.lt and 
s'^ in Eq. (13. 2t upon mq 2uq. (Cf. Figs. 13. 2l and l4. H and discard the Riemann invariants corresponding to 
degenerate gaps). 

Proposition 4.8 After interaction, the boundary between the innermost genus-0 region and the surrouding 
genus-1 regions in cases ( ii) and ( iii) can be written as 

^ = ±(^1 +4ner0 > ^ = ±(^2 > (4-6) 

respectively for the initial portion and for the final portion, where the characteristic speeds are: 4nner ~ '^9^ 
both in case ( ii) and case ( iii); ^-^^g^ = s'^ in case ( ii), and s^^^^^ = s^ in case ( iii). Furthermore, in both 
case (ii) and case (iii) the values of these speeds are: 

4ner = ^0 , = -^O " "0 • (4-7) 

As before, the characteristic speeds can be obtained via degenerate g = 2 calculations (see Appendix IA.4I 
for more details). As before, however, the appropriate values of s'^ , s^ and s'^ can also be obtained from s'^ 
in case (ii) of section |3] upon rescaling go ^ + 5W0 and mq 2uq. More precisely, neglecting degenerate 
gaps, the Riemann invariants that determine the value of ^jj^g^ = s'^ in case (ii) and s^-^J^^^ = Sj in case (iii) 
are obtained from those that determine s'^ in Fig. 13.2b upon rescaling uq —>■ 2uq. (Cf. Figs. 13.2b and l4.1b .c.) 
Similarly, the Riemann invariants that determine the value of ^j^jlg^ = s^ both in case (ii) and in case (iii) 
are obtained from those that determine s'^ in Fig. 13.2b (case (ii) in section |3ll upon rescaling go — > ^0 + ^"o- 
(Again, cf. Figs.E^b andlTTb.c.) 
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Remark 4.9. Equations (I4.5t and (I4.7t also describe the motion of the boundaries of the genus-0 regions 
in the critical case uq = qo, upon proper relabeling of the Riemann invariants. In particular, they predict that 

(2) 

in the critical case the width of the inner genus-0 region is constant (i.e., sl^^^^ = 0). 

The above results, and in particular Eqs. (14.51) and (I4.7t show that the speed of the boundaries between 
regions of genus-0 and genus- 1 after the individual perturbations have come into contact is significantly 
altered as a result of the interaction. In other words, the interaction between different genus- 1 regions results 
in significant bending. This can perhaps be best appreciated in the critical case (i.e., uq = qo, cf. Fig. 14.21 ). 

Finally, we note that, as in the case of one frequency jump, all of the above values for the speeds of the 
boundaries between genus-0 and genus- 1 regions agree very well with the results of numerical simulations 
of the NLS equation, to be discussed in section |6l 



5 Behavior of finite-genus solutions: Arbitrary number of jumps 

Similar calculations to the ones described in section|4]can be repeated for jumps of arbitrary size and/or when 
more than two frequency jumps are present. In particular, from the generalization of the above calculations 
to an arbitrary number of jumps, it is possible to extract some general features. Thus, we consider an initial 
condition with constant-amplitude and with initial frequency expressed in terms of frequency jumps each 
of size Cj and located atx = Xj, j = 1 . . . ,A'^: 

k 

u{x, 0) = -Ctot/2 + ^ Cj for Xk<x<Xk+u (5.1) 
;=i 

where Xq = — oo and X^^i = oo, and where Qot = L^=i Cj = u{oo,0) — m(— oo,0) is the total jump size. The 
constant frequency offset Ctot/2 in u{x,0) does not affect the qualitative behavior of solution, and is chosen 
so that the mean frequency of the solution is zero, i.e., so that the mean position of the solution is constant. 

For simplicity, hereafter we will limit ourselves to consider a symmetric collection of positive frequency 
jumps. That is, we will assume that Cj > V j = I,... ,N and that if a given jump of amplitude C = C* 
exists at X = X^, another jump with the same amplitude exists at x = —X.^,. (Or, in other words, the initial 
condition u(x,0) for the frequency possesses reflection symmetry with respect to the origin.) Since Cj > 
V; = 1 , . . . , A^, the analog of Eqs. d4l1 i holds: 

maxr?(x) = r?(oo) = lim r?(x) = -2qo - Ctot/2 , (5.2a) 

.veR Jf^t" 

minr5(x) = r5(-oo) = hm r5(x) = 2^o + Ctot/2 , (5.2b) 

xeR x^-oo 

Note that again r°(x) and r2(x) are the Riemann invariants before regularization. Owing to Eqs. (15. 2t and 
Definition l4.2l a collection of jumps will be super-critical if Ctot > 2qQ, sub-critical if Ctot < 2qo and critical 
if Qot = 2^0- 

Proposition 5.1 In the case ofN equal frequency jumps each of size C = 2uq, no individual genus-0 regions 
develop if the jumps are individually super-critical, i.e., uq > 2qQ. 

Recall that the size of each of the jumps is 2uo and that the jumps are individually super-critical according 
to Definition 13.11 when the original Riemann invariants overlap at a single location in space, that is, if 
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r2{Xj ) < r^{Xj^) for j = l,...,N. The result applies independently of the number of jumps. More in 
general, however for an arbitrary collection of frequency jumps we have: 

Theorem 5.2 It is possible to obtain arbitrarily large genera for finite times by considering appropriate 
collections of frequency jumps. However, asymptotically as t ^oo^ the only expanding regions in the solution 
are of genus -0, genus-1 and genus-2. 

Proof. The result follows directly from Proposition 12.11 (double sorting property of the characteristic 
speeds). Let us first note that the non-regularized system has two Riemann invariants, and therefore, in- 
dependently of the behavior at intermediate values of x, the regularized system will have only two non- 
degenerate branches as ;c ^ ±00, i.e., /^{x) < r2{x). The only topologically distinct ways to connect these 
two branches are shown in Figs. I5.1h -d (some of these possibilities were also studied in Ref. [10], but not 
for genus-2 cases), where we have neglected the possible presence of one or more shrinking regions such 




Figure 5.1: Qualitative diagrams illustrating topologically distinct cases for the asymptotic values of the Riemann 
invariants after removal of the degenerate portions, together with the corresponding separation between regions of 
genus-0, genus-1 and genus-2: (a, top left) genus-0 region at the center, corresponding to Figs. 13.2b and 14. lb : (b, 
top right) genus-2 region at the center, corresponding to Fig. 14.11 1: (c, bottom left) genus-1 region at the center, 
corresponding to Fig. l3.2h : and (d, bottom right) originated by a super-critical jump surrounded by sub-critical ones. In 
each case, the wavy lines show the branch cuts in the Riemann surface = ny£|^(-Z~ ^ corresponding 

to the local specttum of the Lax operator of the NLS equation. 
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as the outermost genus-2 regions in Fig. l4.1H (see Remark l54l later) . The distinction between the different 
cases in Figs. l5.1h -d is based on the size of the frequency jumps across the pulse. More precisely: 

(i) Fig. 15. Ih is obtained when r^{°°) < r2{—°°), i.e., when the jumps are collectively sub-critical (as in 
Figs.lTIbandlTTt)). 

(ii) Fig. 15. lb is obtained when r°(oo) > r2{—°°) but (x+) < ^^(x^) Vx S M; that is, when the jumps are 
collectively super-critical but every jump is individually sub-critical (as in Fig. l3.2H ). 

(iii) Fig. 15.1b is obtained when there is only one jump, but that jump is super-critical; that is, a single jump 
at X = such that r°(x+) > ^■^(x;^) (as in Fig. l3.2h ). 

(iv) Fig. 15. in is obtained when there is one super-critical jump at the center, surrounded by others jumps 
which are all sub-critical. That is, when there is one jump at x = such that (0+) > r2(0^), sur- 
rounded by other jumps at x = tj such that r^{xj) < r2{xj) Vxy 7^ 0. 

Of course, many variations of these basic configurations are possible depending on the details of the arrange- 
ment of the frequency jumps. In addition, cases (iii) and (iv) (corresponding to Figs. l5.1b .d'). admit a variant 
in which one or more localized "islands" exist between the two outer portions, as shown in Fig. l5.2h .b (cf. 
Fig. 14. Ih ). One of these islands is obtained whenever two additional frequency jumps are inserted which 
are both individually super-critical, that is, when (x+) > r2{x^) for more than one value of x* G M. For 
example. Fig. 15. 2h is produced when there are four super-critical jumps (cf. Proposition I5.lt : Fig. 15.2b is 
produced when there are two super-critical jumps surrounded by two sub-critical ones. From these exam- 
ples it should be clear that the above situations can be easily generalized to cases where an arbitrary number 
of islands is present by adding a proper number of super-critical frequency jumps. Note that the range of 
values for the Riemann invariants covered by each of these islands is disjoint from that of all of the other 
islands. (That is, the islands are not lying side by side. Rather, they are arranged by increasing values of the 
invariants.) Thanks to the double sorting property possessed by the Riemann invariants (namely, Eqs. (12.41) 
in Proposition I2.lt . islands characterized by larger values of the invariants propagate faster to the left, and 
therefore, if several islands exist, different islands will temporarily end up stacked on top of each other 




Figure 5.2: Additional diagrams of asymptotic configurations of Riemann invariants, illustrating the possible presence 
of separate "islands", each of which produces one additional region of genus-2. Figure ls!2h (left): a variant of Fig. lOt 
corresponding four supercritical jumps. Fig. 15.2b (right): a variant of Fig. 15. IB corresponding to two super-critical 
jumps surrounded by two sub-critical ones. 
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upon propagation. It is then clear that by properly choosing the number and size of frequency jumps, it is 
possible to produce regions of arbitrarily high genus for finite values of t. However, the sorting property of 
the Riemann invariants also implies that, eventually, each of the islands will separate from the others. Thus, 
asymptotically in t, at most one island will be present inbetween the top and bottom branches for each value 
of X, implying that corresponding to each value of x there will be at most two gaps. □ 

Remark 5.3. The above decomposition of the solution of the NLS equation (11.11) into regions of genus-0, 
1 and 2 asymptotically with time should be compared to the corresponding result for the Korteweg-de Vries 
(KdV) equation Ut + 6uux + m^xc = 0, for which it was shown that the solution decomposes into expanding 
regions of genus-0 and genus- 1 asymptotically in time.^ 

Remark 5.4. It is important to note that regions of genus-3 may also exist in the solution of NLS for all 
finite values of time. For example. Fig. 15. 3h shows the invariants produced by a configuration with two 
collectively sub-critical jumps surrounded by two individually super-critical ones, while Fig. 15.3b shows 
the invariants produced by a collection of three jumps which are individually sub-critical but collectively 
supercritical. Due to the double sorting property of the Riemann invariants, however, the width of these 
regions decreases monotonically, and tends to a constant asymptotically with time. Therefore, the presence 
of these regions does not invalidate the statement of Theorem 15. 21 Moreover, we expect the width of these 
genus-3 regions to tend asymptotically to zero with time (cf. Conjecture 14. 6t . The situation is similar to 
the case of the KdV equation, where, in addition to expanding regions of genus-0 and genus- 1, shrinking 
regions of genus-2 may also exist asymptotically. 

From the proof of Theorem l5.2l we also have: 

Corollary 5.5 Whenever an overall genus-0 region develops at the center of the pulse as a result of the 
interaction among all frequency jumps, its amplitude is always given by the analogue of Eq. (II. 8I) . namely 
Pmax = j(, (''2(°°) ~ ''i (~°°)) ' similar way as for the single-jump case. 

In other words, the above value is the largest amplitude of a genus-0 region that can be achieved with any 
number N of frequency jumps (cf. Figs. 13.2b and 14. Ib .c). In the case of equal frequency jumps of size 




Figure 5.3: Additional diagrams of asymptotic configurations of Riemann invariantes showing explicitly the presence 
of shrinking genus-3 regions; Fig. l5.3h is the analog of Fig. l5T2b : Fig. 15. 3b the analog of Fig. 15. lb (see text for details). 
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C = 2uo, it is pmax = ?o(l +A^"o/2^o)^- 

Corollary 5.6 Asymptotically as t ^ °°, the solution in the central portion of the pulse will be a region of 
genus-0 iff the jumps are collectively sub-critical, i.e., iff rj(oo) < r2{—°°). 

In the case of equal-size frequency jumps, this condition is satisfied iff mq < 2^o /N . Thus, the maximum am- 
phtude of a stable genus-0 region is always p^ax = 4^0' again as in the single-jump case, and independently 
of the number of frequency jumps. 

All of the results in this work are relative to initial conditions with constant amplitude and piecewise- 
constant frequency. If the initial amplitude and frequency vary continuously, some differences can obviously 
be expected in the solution. Because the quasilinear system of PDEs for the regularized Riemann invariants 
is always strictly hyperbolic, it is stable with respect to small changes in the initial conditions. Thus, even 
though the previous results are relative to the case of discontinuous amplitude and frequency jumps, they 
describe generic features of the solution. Consider for example the dam breaking problem. If the amplitude 
varies continuously from to its maximum value, oscillations appear near the edges of the pulse. (Indeed, 
these oscillations are visible in the numerical simulations.) In the absence of frequency jumps, this scenario 
was studied by Forest and McLaughlin using the Lax-Levermore theory' As long as the transition region 
is narrow, however, the amplitude of these oscillations is small, and therefore the solution will bear a close 
resemblance to the solutions described here. Similar considerations apply if the frequency transitions are 
not discontinuous, as will be discussed in section |6l 

The above results are derived assuming a symmetric collection of positive frequency jumps. The pres- 
ence of negative frequency jumps at the edges or inbetween other jumps would obviously affect the quanti- 
tative details of the picture. We expect that the qualitative features, however, and in particular the asymptotic 
decomposition into regions of genus-0, 1 and 2, to remain valid. 

6 Behavior of finite-genus solutions: Numerical simulations 

Even though the analytical calculations described in the previous sections provide a general picture of the 
behavior of finite genus solutions of the NLS-Whitham equations, there are obviously limits to what can be 
done analytically. In this section and the following one we therefore complement those analytical results by 
presenting numerical simulations of the defocusing NLS equation with small dispersion. 

All the numerical results presented are obtained by numerically integrating the NLS equation (11.11) 
with £2 = 0.1 using a fourth-order Fourier split-step method. We consider an initial condition q{x,Qi) = 
|^(x,0)| exp[j(p(:t,0)/8], where the initial phase ^{x,0) is obtained by integrating Eq. (I5.lt : 



where the function H{x) = \ii x>0 and //(x) = if x < represents the Heaviside unit step. The initial 
pulse amplitude was taken to be a super-Gaussian, namely |^(x,0)| = q^e^'^^/^^^ with M = 40. The 
integration step size and the width AXq of the pulse were chosen to be respectively sufficiently small and 
sufficiently large that none of the numerical results described in this work are affected by them. In order 
to remove the factor £ in front of the time derivative in Eq. (Il.lt . the numerical simulations were set-up in 
terms of the fast scale t' = t /z. Accordingly, all the numerical results will be described in terms of this time 



N 




(6.1) 
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variable. In what follows we describe the behavior of the solutions for varying values of the parameters Xj 
and Cj and for varying number of frequency jumps. We first consider the case of two frequency jumps, to 
numerically identify regions of different genus (Fig. I6.lt and show a few snapshots of the time evolution 
(Fig. I6.2t . We then proceed to describe three sets of simulations: The first set (Fig. 16. 3t shows the effect 
of gradually increasing the size of the frequency jump. The second set (Fig. I6.4t describes the effect of 
increasing the number of jumps. Finally, the third set (Fig. 16. 5t shows the effect of changing the initial 
separation between these jumps. Although numerical experiments were performed with various values of 
amplitude, all the figures in this work except Fig. l6.8h are relative to the case = 1- 

Two frequency jumps. Figure l6Jl shows contour lines of in the (x,f)-plane for different value 

of the jump size Ci = C2 = 2uq =: C. Because of the finite, non-zero value of £, not all of the features 
presented in the previous section are immediately apparent. In other words, because of the finite value of £, 
the solution of the NLS equation is only approximately described by a finite-genus solution in each region. 
Moreover, numerically reconstructing the genus of the solution is a highly nontrivial problem, which is 
outside the scope of this work. Nonetheless, some regions of genus-0 and genus- 1 (whose boundaries have 
been delimited by dashed lines in Fig. I6.lt are recognizable (e.g., regions of genus-0 are often fiat, and 
therefore appear white in the contour plots in Fig. I6.lt . and their location displays a very good qualitative 
and quantitative agreement with the analytical calculations presented in section |4] Note in particular that, 
even though the value of £ is not particularly small, the value of C which corresponds to the critical case is 
exactly that which was predicted analytically, namely C = 2. Figure lOl shows four snapshots illustrating 
the evolution with time of the critical case shown in Fig. 16. lb . 

It is also interesting to note that the hexagonal patterns visible in Figs. l6.1l (in particular, near the center 
;c = in Figs. 16. Ih and 16.1b ') show the interaction of two counter propagating genus- 1 (periodic) waves, 
which produce, as a result, a genus-2 solution of the NLS equation. Note that in the absence of nonlinearity, 
a superposition of two simple waves would produce just a parallelogrammic pattern. Thus, the hexagonal 
pattern is the result of nonlinear phase shifts in the interaction. (Note that the hexagonal pattern can also 
be found in genus-2 solutions of the Kadomtsev-Petviashvili equation describing two-dimensional shallow 
water waves^.) 

Size of the frequency jumps. This set of simulations, some results of which are shown in Fig. 16.31 is 
relative to a two-jump initial condition. We take the two jumps to be equal-size, Ci = C2 = 2uq = C, and 
to be initially positioned at Xi = X2 = ±2.5. Looking at the behavior of the solution as C varies, one can 
clearly observe that the amplitude of the genus-0 portion of the solution increases linearly with increasing 
value of C. More precisely, \q\max = <?o + C/2, in agreement with the analytical calculations presented in 
the previous section (cf. Corollarv l5.5t . See also Fig. l6.7h later. Moreover, in a similar way as in the single- 
jump case, one can see that the maximum amplitude of a stable genus-0 region is |^|max = 2<7o. When uq is 
larger than the threshold value q^, the genus-0 region eventually disappears, as explained in section|4]and as 
described in Fig. 14.2b . Unlike the single-jump case, however, one can obtain larger values of p for limited 
times, as shown in Figs. I6.3b .d. Note how, over these shorter times, the amplitude of the genus-0 region 
is still given by |^|max = ^0 +C/2. Note also, however, that the genus-0 solution closes more rapidly for 
increasingly large values of C, and, above the threshold C = 4 (i.e., C = ^0) does not open at all, at least for 
a fixed value of AX. These results are in very good agreement with the analytical calculations described in 
the previous sections (e.g., see Proposition 15. H and Corollaries 15.51 and 15 . 6t . 
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Number of frequency jumps. In this set of simulations, some of which are shown in Fig. 16.41 the number 
N of frequency jumps was varied, while the size of each jump and their initial spatial separation are kept 
fixed respectivaly at C = 2 and AX = 5. From Fig. 16.41 it is apparent that the amplitude of the genus-0 part 
also increases linearly with the number of jumps. More precisely, |^|max = Qo +NC/4, in agreement with 
the analytical calculations presented in the previous sections (cf. Corollary 15.51) . See also Fig. 16.7b later. 
The genus-0 region, however, will eventually close up if the to to/ jump exceeds the same threshold as in the 
previous case: Qot = NC = 4qo (cf. Corollarv l5.6t . For increasingly larger numbers of jumps, the amplitude 
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Figure 6.1: Contour plots of numerical simulations of the NLS equation ( II. H for a two-jump initial condition and for 
different values of the jump size parameter C: (a, top left) C — 5, corresponding to Fig. 14. lb and Fig. 14.2b . showing 
two jumps that are both individually super-critical; (b, top right) C ~3 corresponding to Figs. l4.li :.d and Fig. 14.2b . 
showing two jumps which are individually sub-critical but super-critical in combination; (c, bottom left) the case C — 2, 
corresponding to Fig. 14.2b in which the two jumps are collectively critical; (d, bottom right) C = 1, corresponding to 
Fig. 14. lb and Fig. l4.2H and showing a case in which the jumps are collectively sub-critical. In each case, the dashed 
lines demarcate the recognizable regions of different genus. In particular, the speeds of the boundaries between genus- 
and genus- 1 regions are in very good agreement with the analytical results described in the previous sections. 
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of the genus-0 region follows the same law as above, namely |^|max = +NC/A, and can temporarily 
achieve very large values. The genus-0 region also tends to close up more rapidly, however, and, above a 
certain number of jumps, does not open at all (at least for a fixed value of AX). 

Spatial separation between frequency jumps. In this set of simulations we look at the effect of increas- 
ing the spatial separation AX between the frequency jumps. Figure 16.51 is relative to a two-jump initial 
condition with C = 2. The consequences of changes in AX are not as obvious as those of changes in the 
number of jumps or the jump size. One can see, however, that increasing AX has two main effects: on 
one hand, it lenghtens the time scales over which the evolution occurs, since the various regions need to 
travel for longer distances to come into contact and to interact; at the same time, it determines the width 
of the stable genus-0 region in the critical case. In fact, numerical results show a remarkable fact, namely 
that the width Xcx. of the stable genus-0 region depends linearly on the initial spatial separation AX between 
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Figure 6.2: Numerical simulations of the NLS equation illustrating the evolution with time of a critical two-jump initial 
condition. The frequency jumps, located at X2 = — Xi =5/2 have amplitude Ci =€2= 2, so that their interaction 
produces a genus-0 region of constant width in the center (as shown in Fig. l6.lt ). The solution is shown at: (a, top 
left) /' ~ 2, (b, top right) t' = 3.5, (c, bottom left) t' = 5, (d, bottom right) t' — 10. This case corresponds to the contour 
plot in Fig. l6.lt and the bifurcation diagram in Fig. l4.2t . 
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Figure 6.3: Numerical simulations of a two-jump initial condition with ziX = 5, for different values of the jump size 
Ci = C2 = C: (a, left) C = 1; (b, right) C 7/2. The critical case C = 2 is also shown in Fig.lOll. In Fig.lQi 
the solution is subcritical (the genus-0 region expands forever, see Fig. I4.2t l). and the solution is shown at t' ~ 10, 
whereas in Fig. 16.3b the solution is supercritical (see Fig. 14.2b ). and it is shown at the time when the genus-0 portion 
has maximum width, t' ~ 5.2. The horizontal dashed line at \q\ ~ 1 represents the amplitude of the initial condition, 
and the other horizontal dashed line identifies the amplitude of the genus-0 region. Note that: (i) the amplitude of the 
genus-0 portion of the solution increases linearly with increasing value of C. (ii) for increasingly large values of C, the 
genus-0 solution tends to close more rapidly, and, above a certain threshold, does not open at all. 
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Figure 6.4: Numerical simulations for different numbers of frequency jumps, with the size of each frequency jump 
fixed at C, — 2, i — 1, . . . ,A^: (a, left) N ^ 3; (b, right) N = 5. In both cases the solution is shown at t' = 10. See 
also Fig. 13. lb . illustrating the case = 1, and Fig. l6.2t l. illustrating the (critical) case N = 2. Note that for = 1 the 
central genus-0 region expands and for = 2 it has a steady width, whereas for A^ > 2 it exists only for limited times. 
In all cases the amplitude of this central genus-0 portion increases linearly with the number of jumps. However, for 
larger number of jumps, this genus-0 portion tends to close more rapidly, and above a certain number A^ it does not 
open at all, at least for a fixed value of AX. 
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the jumps. In the case of two frequency jumps, one simply has Xoo = AX /2. (Cf. Figs. 16.1b . 16.21 and 16.51 ) 
A similar relation however holds independently of the number of jumps. For a critical arrangement of an 
arbitrary number N of equal-size, equally spaced positive jumps, it is Xoo = {N — 1) AX/2, where AX is the 
initial separation between consecutive jumps. (Cf. Fig. 16.61 and note that {N — I) AX is just the distance 
among the farthest two frequency jumps.) These relations are verified in Figs. l6.8h .b. 

Summary. Let us recapitulate for convenience the main results of the numerical experiments discussed in 
the above paragraphs. Given N positive frequency jumps with individual sizes Cj and sepai^ations AXj, we 
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Figure 6.5: Numerical simulations for different separations between jumps, here with a critical arrangement of two 
frequency jumps of amplitude Ci = C2 = 2: (a, left) AX = 2, (b, right) AX = 16. See also Fig. l6.2t l. which shows the 
case AX = 5. In Fig. l6.5h the solution is shown at t' = 10, in Fig. 16.5b at f' = 32. Increasing the value of AX has two 
main effects: it increases the maximum width of the genus-0 portion, and it simultaneously lengthens the time scales 
over which the evolution occurs. See also Fig. l4.2t . 




Figure 6.6: Critical arrangements of more than two jumps: (a, left) N ^3 and C = 4/3; (b, right) N — 4 and C = 1. In 
both cases qo = 1 and AX = 4. 
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can observe that: 

(i) The value of the asymptotic genus in the central portion of the pulse is only determined by the total 

jump Ctot = Ci H VCn (or, equivalendy, by the size per jump C = Ctot/N in the case of A'^ jumps of 

equal size). More precisely, the solution will be genus-0 iff Qot < ^o, in very good agreement with the 
analytical results (see Corollarv l5.6t . Of course, the detailed behavior of the solution for finite times 
depends on the precise values of all the Cj and AXj. 

(ii) If a genus-0 region is desired which expands forever at the center of the pulse, its height is limited 
by the critical value of the cumulative jump Ctot- More precisely, |^|max = 2^o, independently of the 
number of jumps. For equal-size jumps, this value is obtained for C = Aq^/N, in very good agreement 
with the analytical results (see Corollarv l5.5l) . 

(iii) It is possible to achieve genus-0 regions of much larger amplitudes over limited times by employing a 
sufficiently large number of jumps. More precisely, |^|max = ^0 + ^01/4 (i.e., |^|max = ^0 +A'^C/4 for 
equal-size jumps), as shown in Fig. 16.71 and again in very good agreement with the analytical results 
(see Corollarv l5.5t . 
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Figure 6.7: (a, left) The numerically calculated amplitude of the genus-0 region produced by an arrangement of 
equal-size frequency jumps as a function of the indidual jump size C. Triangles: = 1; squares: N ^ 2; circles: 
N ~3. For comparison, the solid, dashed and dot-dashed lines show the straight lines y = {N — \ )C/A. (b, right) The 
amplitude of the genus-0 region produced by an arrangement of frequency jumps each of size C = 2, as a function of 
the number of jumps. Squares: numerical simulations; dashed line: the \mty— l+N/2. 
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Figure 6.8: Width of the genus-0 region resulting from the interaction of a critical set of frequency jumps as a function 
of the spatial separation between them: (a, left) N = 2; solid line: the straight line y = AX/2; squares: simulations 
with ^0=1 and C = 2; triangles: simulations with qo = 2 and C = 4; stars: simulations with = 3 and C = 6. (b, 
right) N = 3 (squares) and N = 4 (triangles), with q'o = 1 in both cases; solid line: the straight line y = {N — l)AX/2. 
Note that, since the characteristic period! of the rapid oscillations is (9(e), it is only possible to measure the width of 
the genus-0 region with precision Ax ~ ±e. 
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(iv) For a critical arrangement of an arbitrary number of equal-size, equally spaced positive frequency 
jumps, the width of the stable genus-0 region is half of the initial separation between the farthest two 
frequency jumps. This fact has not been demonstrated analytically, but is nonetheless very solidly 
supported by the numerical results, as shown in Fig. 16.81 

In the next section we discuss how these results (and in particular items iii and iv) can be used to generate 
intense, ultra-short optical pulses. 

7 Applications: Generation of intense short optical pulses 

The analytical and numerical results described in the previous sections can be used to design an initial 
condition that results in a genus-0 region of arbitrarily high amplitude at any desired distance down the fiber 
by employing a sufficiently long continuous wave (CW) pulse at the outset, as we now describe. (Recall that, 
for fiber optics x is the retarded time and t is the propagation distance along the fiber.) Of course, the same 
techniques can also be applied in any physical context where the dispersionless limit of the NLS equation 
is relevant. The basic building block consists of two frequency jumps of critical size, which generate a 
genus-0 region of constant width, as shown in Fig. 16.1b . By including additional frequency jumps to both 
sides of this basic pair, one can then temporarily produce an overall genus-0 region of higher amplitude. 
This basic process is illustrated in Fig. 17. II which shows the time evolution of a four-jump initial condition. 
A contour plot of this solution is shown in Fig. 17. 2h . We now discuss some practical issues related to the 
generation of these high-intensity pulses: namely, how to control their amplitude, the distance along the fiber 
at which they are produced, their temporal width, and post-processing through filtering in order to eliminate 
the high-frequency oscillations. 

Amplitude, location and width. By appropriately choosing the number, amplitude and location of the 
additional frequency jumps, it is possible to precisely specify the amplitude of the overall genus-0 region 
as well as to prescribe its width and the distance along the fiber at which it is generated. This is illustrated 
in Fig. 17.21 which contains contours of four different initial conditions. Figure 17. 2h corresponds to the 
four-jump solution shown in Fig. 17. 11 the two central frequency jumps are critical (C2 = C3 = 2), while the 
single, sub-critical frequency jumps on either side of them {C\ = C4 = 3) provide two expanding genus-0 
regions which go to interact with the main one resulting in an overall genus-0 region of amplitude |g|max = 
^0 + Ctot/4 = 7/2, once more in perfect agreement with the analytical prediction. 

With similar arrangements of four jumps, one can obtain any value of amplitude |g|max < 4, at which 
point the two side jumps become critical, and do not generate a genus-0 region anymore. This limitation 
on the pulse amplitude can be overcome, however, by considering a larger number of frequency jumps. 
Figures I7.2b .c are both relative to a six-jump initial condition with all jumps having amplitude Cy = 2, so 
that each pair of adjacent jumps generates a genus-0 region of constant width. These regions then interact 
to create an overall genus-0 region of amplitude 4. In Fig. 17.2b . the jump locations — X3 = X4 = 2, —X2 = 
X=, = 15, —X\ = Xe = 20 were chosen so that the overall genus-0 region arises near t' = 15. In Fig. 17.2b . 
instead, the side jumps —X2 = X5 = 20, — Xi = Xg = 35 were positioned further away, so that the genus-0 
region emerges near t' = 25 and is maintained for a longer interval. Finally, Fig. 17.21 1 shows an eight-jump 
initial condition, with equal-size jumps Cj = 2, generating an overall amplitude |(7|max = 5. Once more, 
the jump locations were chosen so that this region emerges near t' = 25. A convenient way to choose the 
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jump locations is to remember that each frequency jump moves with velocity corresponding to the average 
frequency across the jump. 

Finally, we note that the width of the central genus-0 region can be easily adjusted by changing the 
temporal separation between the two critical frequency jumps, by virtue of the linear relation between these 
two quantities which was discussed in section |6l 

Smooth frequency transitions. The remarks at the end of section |5] about the general stability of the 
analytical results with respect to small changes in the initial conditions also apply with regard to the dis- 
continuity in the frequency jumps. That is, we expect that the general picture described here will be an 
approximate description of the pulse behavior even if the frequency transitions were continuous, as would 
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Figure 7.1: Numerical simulations of the NLS equation with a supercritical four-jump initial condition. The frequency 
jumps, located at X4 = —Xi = 24 and X3 = —X2 = 2, have amplitude C2 = C3 = 2 (so that the interaction of the two 
center jump produces a critical genus-0 region in the center) and Ci —C4 — 3, so that the two individual sub-critical 
jumps open up a genus-0 region which moves towards the center of the pulse. The solution is shown at: (a, top 
left) t' = 2, (b, top right) t' = 6, (c, bottom left) f' = 16, (d, bottom right) t' = 24. A contour plot of this solution is 
shown in Fig. l7.2h . 
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be the case in any experimental setup. To test this prediction, we performed numerical simulations in which, 
instead of Eq. (I5.lt . the initial frequency was chosen according to 

1 ^ 

u{x,0) = -Y,Cjtanh{x-Xj), (7.1) 

which implies (?{x,0) = {I /2)'£^^i Cj log {cosh{x-Xj)) . The comparison of the numerical results for 
discontinuous and smooth frequency transitions shows remarkable agreement between the two cases, with 
regard to both the bifurcation diagrams and the the amplitude of the genus-0 regions. As an example, 
in Fig. 17.31 we show numerical simulations of two frequency jumps with the same parameter values as in 
Figs. l6. Ib -d. Note also that the hexagonal pattern near the center ;c = in Fig. l7.3h shows a genus-2 solution 
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Figure 7.2: Contour plots of the numerical solution of the NLS equation with various choices of frequency jumps: (a, 
top left) four jumps, with parameters corresponding to the solution shown in Figs. I7.1h -d: (b, top right) six jumps, 
positioned at — = = 2, —X2 = = 15, —Xi =X(, = 20 and all with amplitude C — 2; (c, bottom left) six jumps, 
with same parameters as before except —X2 = X5 = 20 and —Xi = X(, = 35; (d, bottom right) eight jumps, each with 
amplitude C = 2 and positioned at -X4 = X5 = 4, -X3 = 20, -X2 X7 = 35, -Xi = 50. 
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of the NLS equation (see also Figs. 16. II) . The evolution of the Riemann invariants corresponding to Fig. l7.3h 
is given by a smooth- version of Fig. lSTTI 

Filtering. The high-frequency oscillations that appear on each side of the genus-0 region may be undesir- 
able for practical purposes. These oscillations can be effectively removed by appropriately filtering the pulse 
at the fiber output, resulting in a smooth, high-amplitude short optical pulse. For example, post-processing 
with a Lorentzian filter with frequency response F((o) = {-s/R/n) / (1 +/?(0^) with /? <C 1 effectively "wipes 
out" the genus- 1 and genus-2 oscillations. A drawback of this method however is that, in the case of a 
narrow pulse, the conversion efficiency (i.e., the amount of energy in the initial CW pulse that is preserved 
in the final high-amplitude pulse) is rather low after filtering. A higher conversion efficiency can be obtained 
either by employing an Erbium-doped fiber amplifier or by increasing the temporal separation between the 
two central frequency jumps, so as to obtain a wider genus-0 region. 

Parameter scaling. As a concluding remark we note that, even though almost all the figures in this work 
describe solutions produced by initial conditions with amplitude qo= 1, all of these numerical results can 
easily be related to different values of initial amplitude thanks to the scaling symmetry of the NLS equation. 
Namely, ifq{x,t) —>■ aq{x,t), the results in this work still apply upon rescaling x —>■ x' = x/a, t —>■ t' = t /a-^ , 
implying p{x,t) —>■ p'{x' ,t') = a^p{ax' ,a^t') and u{t,x) — > u'{x' ,t') = au{ax' ,a^t'). Importantly, the same 
scalings can also be used to relate the results to choices of spatial or temporal units different from those 
discussed in Appendix lA.il 
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Figure 7.3: Numerical simulations of NLS with smooth frequency transitions, with two frequency jumps located at 
X2 — —Xi = 2: (a, left) C ~3; (b, center) C — 2, (c, right) C = 1. These figures should be compared to Figs. l6.1b -d. 
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Appendix 



A.l Nondimensionalizations and scalings 

It has been known for more than thirty years that the propagation of a coherent Ught pulse in optical fibers 
is governed by the NLS equation, which, in physical units, is^H^ 

BE 1 ..d'^E , 

where E{t,z) is the slowly varying amplitude of the complex envelope of the electric field of the pulse, 
t = tiab ~z/cg is the retarded time and z is the propagation distance. Here, Cg = l/k' is the group velocity, 
^(co) the propagation constant, k" = (f-k/dap- the dispersion coefficient and y the nonlinear coefficient. In 
the following, we will use the typical value 7= 3W^^km^^ 

Introducing the dimensionless variables q = E j \fPl, T = t /t.^ and Z = z/Zn\, where and are some 
characteristic power and temporal duration, and where z* = Zni = 1 / (Y^*) is the nonlinear length (the char- 
acteristic distance over which nonlinear effects take place), Eq. dA.ll ) becomes 

The dimensionless dispersion coefficient d = k" /k'l = Zni/zdisp (with k'l = tl/zn\) quantifies the relative 
importance of nonlinear and dispersive effects in Eq. dA.ll) . In this work we are interested in particular to 
situations where < <i ^ 1. This can be achieved in two different ways: constant dispersion or dispersion 
management. Note that the dimensionless dispersion coefficient d in Eq. (IA.2t is obtained equivalently as 
d = k" /k^ or as d = D/D^, withD = —{2%c/}?)k", with = —{2%c/}?)k'l, and where D is the dispersion 
coefficient in ps/(nm-km). Note llicj)^ = 0.78 (nm-ps)^^ at at X = 1.55/7m. 

Typical units. We first consider the case of fibers with constant dispersion. A unit power = 2 W (which 
can be obtained using commercially available high-power Erbium-Doped-Fiber-Amplifier lasers) implies 
Zni = 0.16 km. Then a characteristic time unit of = 50 fs (which is appropriate in order to describe short 
optical pulses) yields k'l = 1.5ps^/km. Then the dimensionless dispersion coefficient d in Eq. (IA.2t is d = 
D/D^, withD* = — 1.2ps/(nm-km). A value of (/= 0. 1 (such as the one used in all the simulations presented 
in this work) can then be obtained by employing a fiber with dispersion coefficient D = — 0.12ps/(nm-km). 
Alternatively, the same effects can be observed for lower-amplitude, longer pulses over longer propagation 
distances: for example, with a unit power = 20 mW, and a unit time = 0.5 ps, the unit distance would 
be z* = 16 km. 

Fibers with a larger nonlinear coefficient (such as a dispersion-shifted fiber or a photonic-crystal fiber, 
for example) would have a correspondingly smaller nonlinear length Znh which would make it possible 
to observe the same nonlinear effects over shorter distances. Also, for comparable values of and P^, a 
higher nonlinearity would result in a larger value of k'l, which means that fibers with correspondingly larger 
dispersion coefficients may be used. Alternatively, a higher nonlinearity would make it possible to achieve 
the same value of kl (and hence d) with shorter temporal scales and/or with smaller peak powers P^ . 

It should be noted that Eq. dA.ll) neglects the effect of fiber loss, and therefore the peak powers listed 
above should be taken as the average of the pulse powers over the whole transmission span. This is not 
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a serious issue if the total propagation distance is much shorter than the characteristic distances for fiber 
absorption, as with the sample parameters chosen earlier (for standard telecommunication fibers, a power 
loss coefficient of 0.2 dB/km is typical), or if the effects of damping can be minimized by placing one or 
more (Erbium-doped and/or Raman) fiber amplifiers inside the transmission span, since it is well-known 
that the average dynamics in systems with loss and periodic amplification are still governed by the NLS 
Eq. (Il.lt with constant coefficients.^ We should also note that quasi-lossless propagation of optical pulses 
over several hundred kilometers has recently been achieved by making use of two bidirectional Raman 
pumps plus fiber Bragg gratings.''^' Finally, note that Eq. (IA.lt also neglects higher-order dispersion. As 
shown in Ref. l25l . third-order dispersion can affect pulse behavior. In order for the results described in this 
work to apply, therefore, it is necessary that the fiber dispersion be approximately constant throughout the 
bandwidth of the pulse, which means that third-order dispersion should be small. 

Dispersion management. If realizing small values of <i in a stable way should require an unpractical 
level of control over fiber dispersion, the problem might be circumvented by making use of dispersion 
management. With dispersion management, a very small value of £ can be achieved as an average between 
opposite the dispersion coefficients of fibers with alternating signs of dispersion, even if the local values of 
dispersion are not individually small. 

The propagation of optical pulses in systems with dispersion management is still described by the NLS 
Eq. (IA.2I) . where however the dimensionless dispersion coefficient d in Eq. (IA.2I) is replaced by a periodic 
function d{z). The specific choice of d{z) is called a dispersion map. The relative effect of the periodic 
dispersion variations can be quantified by the reduced map strength parameter s which is defined as' ' ' % = 
{\d{-) — d\) Za, where now d = {d{ ■ )) and the average is taken over the period of the dispersion map, Za- 
As long as the strength s of the dispersion map is small, it is possible to average over the rapid variations of 
the pulse profile originating in the NLS equation with periodic coefficients. Indeed, it is well-known that, 
just hke with loss/amplification, the NLS Eq. (IA.2I) describes the evolution of the leading-order portion of 
the pulse envelope'^ in a system with a moderate amount of dispersion management. That is, the leading 
order part of the pulse satisfies again the NLS equation (IA.2t . except that the parameter d now represents 
the average dispersion. Therefore, the use of dispersion management might be desirable because it allows 
one to more easily obtain small values of effective dispersion than for a system with constant dispersion. 
Because on average the system is still governed by the NLS equation, the results presented in this work will 
be preserved as long as the dispersion map strength s is small. If the map strength s becomes large, however, 
the leading order dynamics are described not by the NLS, but rather by a nonlocal evolution equation of 
NLS-type'^ Determining the behavior of the solutions with small J for this type of system is a highly 
non-trivial task, which is beyond the scope of this work. 

A.2 Genus- 1 solutions of NLS and Whitham's averaging method 

Here we briefly review the construction of genus- 1 solutions of the NLS equation and the derivation of the 
genus- 1 NLS-Whitham equations, with the aim of clarifying the connection between the Riemann invariants 
and the spectrum of the Lax operator. For a detailed construction of higher-genus solutions involving theta 
functions, see Refs. 151 (161 . 

In order to construct traveling wave solutions of the NLS equation dl.ll) . it is convenient to introduce the 
fast space and time scales x' = x/z and t' = t /z, and write Eq. dl.ll) as iq^ — ^qx'x' + = 0- We then look 
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for solutions of the form 

=e'[^^+(^-'+^)''](a + /(y + c/)) , (A.3) 

with a,c,X £M. and /( • ) a real function of its argument. Substituting Eq. (IA.3I) into NLS and integrating 
once, one obtains 

(/)'=/- cJi/ + ^if - 03/ + 04 =: Kf) , (A.4) 
where the prime denotes differentiation with respect to x', Gi = —a, O2 = 4(3a^ — X + c^/2), 03 = —aOj, 
O4 = C and C is the integration constant. It is convenient to represent the right-hand side of Eq. (IA.4I) as 

Kf) = (/ - n ) (/ - r2) (/ - r3) (/ - r,) , (A.5) 

where the roots ri , . . . , r4 are related to the coefficients in Eq. (IA.4I) via the relations 

aj = n +r2 + r3 + r4, 02 = n?'2 + HO + n?'4 + ?'2?'3 + ''2''4 + '■3''4 , (A.6a) 
03 = r\r2ri, + rir2r4 + riri,r4 + r2r3r4 , 04 = rir2ri,r4 . (A.6b) 

Integrating formally Eq. iAAl one has J df/ \/Ji{f) = ±(/ ~-^o)- integral then defines the solution 
of NLS in terms of elliptic functions. Note however that real solutions are confined to those values of / 
such that ^(/) > 0. It is then clear that the values of the four roots ri , . . . , r4 (which will be the Riemann 
invariants when performing the Whitham averaging described in the next subsection) determine the edges 
of the forbidden regions, or gaps, which are the values of / for which < 0. Several situations can occur 
depending on the values of the parameters in Eq. (IA.4I) : it is relatively easy however to see that bounded 
solutions only exist when all four roots of ij{f) are real. In this case, labeling the roots as ri < r2 < ri, < r4 
without loss of generality, f{x') is given by the elliptic sine: f{x') = k sn(x',/:), where the elliptic modulus 
k G [0, 1] is given by k^ = {r^ — rT,){r2 — r\)/[{r4 — r2){rs — ri)]. The elliptic parameter'^ is m = k^. The 
value of / oscillates between r2 and r^, and the spatial period of the solution is given by 4K{k), where K{k) 
is the elliptic integral of the first kind.'^' In the limiting case r4 = r3 = —ri = —r2, the period of the solution 
tends to infinity, and one obtains the dark soliton solutions of NLS : for example, when ^ = it is a = 
and f{x') = X t&r\h(kx'). Note that that the solution of NLS is completely specified up to an overall phase 
in terms of the four roots ri , . . . , r4. For a more detailed discussion of single-phase solutions of the NLS 
equation in the focusing as well as defocusing regimes, as well as their stability, see Refs. I8l l24l . 

It is important to note that, for any finite value of the small parameter £, the presence of 8 in Eq. (Il.lt can 
be taken into account by simply replacing {x,t) by the fast variables (x' = x/£,t' = t/z). This phenomenon 
is the origin of the high-frequency oscillations in the solution of NLS with small dispersion coefficient, 
discussed throughout this work. The nontrivial issue, however, is to understand the qualitative behavior of 
the solution in the weak dispersion limit, as formulated in section ^ 

Let us then briefly outline how to obtain the genus-1 NLS-Whitham equations. It is well-known that the 
NLS equation admits an infinite number of conservation laws: 

dF„ dG„ 



dt' dx' ' 

n = I, ... ,00. For example, the first four conserved densities are 



(A.7) 



Fi = \q\\ (A.8a) 



F2 = i{q*q^ - qq^,) , (A.8b) 
F^ = \q,>\^ + \q\\ (A.8c) 
F4 = i{q*q:ex'x' - qql'x'x') - 3i1^P('?*'?.r' - qq*x') ■ (A.Sd) 
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Then we consider small modulations of the genus- 1 solution given by the elliptic function, in the sense that 
the deformation due to the modulation is of order £ over one period of the wave. Following the standard 
perturbation method of multiple scales, we introduce the slow time and space scales x = zx' and t = zt' . We 
now have two sets of scales, consisting of a set of fast scales for the oscillations and a set of slow scales for 
the modulations. With those variables, we have 

3 3 3 3 3 3 
3^^3^^^3^' df^dt'^^di' 

Then, expanding the conserved quantities as F„ = F° + zF^ -\ and G„ = G° + zG\ + ■■■ and taking 

the average of the conservation laws (IA.7t over the fast scale x', at order £ we obtain the NLS-Whitham 
equations in conservation form: 

mijj^_ (A.9) 

dt dx 

Here (F^) and (G^) are the average of the leading order conserved densities and fluxes associated with the 
solution of the NLS equation, and the average of a periodic function f{x!,t') = /(6) is given by 

. L/2 1 ^ 

(/) = 7 / f{x,t')dx' = jmdQ, (A.IO) 

^ -L/2 

where L is the spatial period and 6 = (x' + ct')/L. For the genus-1 solutions described above, L = AK{k). 
Note that the averaged quantities and {G^) are functions of the spectral parameters r\,... ,r^, Thus 
the modulations are expressed in terms of the slow motion of the spectral parameters r\,. .. ,r^. This also 
means that we need the first four equations from ( IA.9t . By computing the averages, one finds that the 
variables r\,. .. are precisely the Riemann invariants for the system. Then, by writing Eqs. ( IA.9t in 
diagonal form, one finally obtains the genus-1 NLS-Whitham equations (12. 3t with g = \. (That Eqs. (IA.9t 
are diagonalizabile is not obvious, but always true for all values of genus; see section lA3l f or further details.) 



A.3 The NLS-Whitham equations 

Here we briefly review some well-known results for the NLS-Whitham equations. For more details on the 
finite genus solutions of the NLS equation, the averaging process and the NLS-Whitham equations, we refer 
the reader to Refs. |5l|25l and references therein. 



The genus-g NLS-Whitham equations can be obtained by extending the approach described in sec- 
tion to hyper-elliptic solutions of NLS and 2g + 2 conservation laws It is more convenient, how- 
ever, to take advantage of the integrable structure of the NLS equation, as described in Ref. [25 1. Recall that 
genus-g solutions of NLS are associated to the Riemann surface R:w^ = i^g{z), where iJtg{z) = n^lt^(z — 
(cf. Eq. d2!2l ). It was shown in Refs. L12||l3l|22l|29 1 that the NLS-Whitham equations (IX9l can be eleganfly 
written as 

3(0i 3(02 

where COi and 0)2 are the meromorphic (Abelian) differentials of the second kind associated with the Rie- 



38 



mann surface R: 



(Oi 



0)2 = 2 



_ 1 
" 2 

z + 



1 + 



t/z= (l + 0(l/z))c/z, 



(A. 12a) 



z^+2 - iaiz^+' + i (a2 - ^a?)z« + yiz«- ^ + • • • + Y,-iz + Yg 



t/z= (z + 0(l/z))t/z, 

(A. 12b) 



where Oi,. .. ,03 are the generalization of Eqs. (IA.6I) . 
01 = 



I 

i<;<2g+2 



l<j<k<2g+2 



03= 52 

l<!<7<<:<2g+2 



(A. 13) 



(03 will appear later on) and where the coefficients a^- and are given by the solution of two g x g linear 
systems of equations 



1 



ai/f ^ + ^ (02 - ^0?)/J + Yi/f ' + • • • + Yg-i/j + Yg/° = 0, 



(A. 14a) 
(A. 14b) 



j = 1, . . . ,g, with the real- valued functions given by 



1 nj+2 Ji 

/ 

2 r2,+i V-Afglz) 



<iz. 



(A. 15) 



Equations (IA.14t express the requirement that the A-periods of (Oi, (O2 be zero: that is, COi = o>2 = 0, 
where the aj are the A-cycles on the Riemann surface R. 

Evaluating the residue of 0)1,2 at z = in Eq. (lA.llt . one obtains the genus-g NLS-Whitham equations 
in Riemann invariant form: 

-=.,(n,...,r2,+2)3^, 

= 1, . . . ,2g + 2. The characteristic speeds are given by i'/t(?'i, . . . ,r2g+2) = sijk), with 



(A. 16) 



rs+^ - i 0irg+' + l(02-^0f)rg+Yi?-g-i + ---+Yg-ir + Yg 
rS+i - \GirS + a\rS-'^ ^ httg-ir + a^ 



(A. 17) 



Since the NLS-Whitham equations (IA.16I) describe the evolution of the Riemann invariants ?"i , . . . , r2g+2 
with respect to the slow time and space variables, and since the Riemann invariants are also the branch 
points of the Riemann surface corresponding to finite-genus solutions of the NLS equation, it follows that 
the NLS-Whitham equations describe the slow modulation of finite-genus solutions of the NLS equation. 

Note that the abelian differentials in Eqs. (IA.12I) can also be expanded in terms of the conservation laws 
of the NLS equation: 



(Oi 



n=\ 



+ 1 



(O2 



^+1 



dz, 



(A. 18) 
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where and {G^) are the conserved densities and fluxes of the NLS equation, averaged over the fast 
phases. The average of a quasi-periodic function of ^ phases f{x',t') =/(6i , . . . ,6g), with 6yt = {x' + Ckt')/Ljc 
and Lyt being the period for the phase Q^, is given by 

(/) = lim J /' fix'/) dx' = }-- • 1/(61 ,...,Qg)dQi---dQg, 

L^oo L _^i2 

where we have used the ergodic assumption, i.e., that the g phases produce a dense orbit on the Jacobian 
variety of the Riemann surfacel^^. In this framework Eqs. (lA.l It then yield Eqs. ( IA.9t for « = 1 , . . . , 2g + 2. 
For example, the first few densities and fluxes are 

(F0) = (p) = -ia2 + T^a? + iai, 

(/f) = (pw) = -|a3- 10102 + ^C)i + |aiai + ia2, 

(G?) = (F") = -^03 - 1^102 + ^a? + iyi , 

{Gi) = (P"^ + 5P^> = + iCJias + ^cj^ - l^o\o2 + + |aiYi + iyj , 

and Eq. ( 11.51 ) is obtained from Eqs. ( IA.9I ) when g = 0. 

Upon evaluating the integrals in Eq. (lA. 151) and taking the appropriate limits as 8 — > 0+, Eq. (IA.17t can 
be used to calculate the value of all the speeds listed in sections 15] and |3] as we show in the next subsection. 
It should be noted that the general solution of the NLS-Whitham equations (lA. 16t can be written in terms 
of the Hodograph transformation 

X - st{ri r2g+2)t = Wkin,..., rjg+i) , (A.20) 

k= 1 , . . . , 2g + 2, where the functions are expressed in terms of the initial data via hyperelliptic integrals. 
When inverting Eqs. (IA.20t to find the solution q{x,t), the locations where the Jacobian of the transformation 
vanishes determine the boundary between regions of different genus^lEl, xhis information was used in 
Ref. 1 31 1 to determine the location between regions of genus-0 and genus- 1 in some specific situations, and 
the method was used in Ref. 1 17 1 to determine the generic evolution of an arbitrary initial datum for the 
Korteweg-de Vries equation. It should be noted, however, that the method cannot be easily applied for the 
situations considered in the present work, since the transformation x q{x,0) is not invertible in the case 
of piecewise constant initial data, which are the ones of interest here. 

A.4 Calculation of some characteristic speeds 

We now briefly describe some calculations regarding the characteristic speeds of the regularized Riemann 
invariants for the NLS-Whitham equations. The results of these calculations and similar others were summa- 
rized in sections |3]and|4] Since the calculations are rather tedious however, we will limit ourselves to present 
two examples. For more details on similar types of calculations we also refer the reader to Refs. |7 25 1. 
When performing these calculations, it is useful to note that the characteristic speed in Eq. (IA.17t is the 
ratio of two meromorphic differentials which becomes the ratio of two polynomials (the denominator of the 
differentials). In the case of step initial data, some common factors then arise, and the resulting cancellations 
produce a finite result. 
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Genus-1 calculations: Equations O.lt and i4.5\ . Our first aim is to calculate in case (i) of section |3l 
as shown in Fig. I3.2h . This is done by evaluating ST,{ri,. . . ,r4) = ^(rs), with s{r) given by Eq. (IA.17I) 
with r4 = —r\ = uq + 2qo, r2 = — mq + 2go and rj, = r4 — £, and then taking the limit £ 0+. In general, 
Eqs. (IA.14t yield the coefficients tti and yi for the g = I case as 

/?-^ai// 7i^-iai/f + l(a2-H)// 

«i- To > Yi- To • 

^1 ^1 

In our case, Eqs. (lA. 13t also give a\ = r2 + r4 = 4^o + £ and 02 = r2r4 — r^ = — 2mo(mo + 2^o) +£("0 — 2^o)- 
In general, upon making the substitution z = rs + A sin^ 6 with A = r4 — r3. the integrals /{',.•• 5^1^ take the 
convenient form: 

(r3 — ri + A sin^ 0) {tt, — r2 — A sin^ 6) 



In our case, expanding the power in the numerator and substituting the explicit value of the invariants, the 
integrals can also be written as: 

k / h\ It/2 ein2'"fl 

= I («o + 2^0)'^-'" e'" / = . (A.23) 



m=0\^J ' ^(2(Mo + 2<7o)-£COs2e)(2Mo-£COs2e) 

Then, evaluating the integrals in Eq. (IA.23t . substituting the results and Eq. (IA.21I) into Eq. (IA.17t and 
taking the limit £ — > 0+, one then obtains the desired result, namely Eq. (I3.lt . Note that, even though 
limE^o+ ^1* = the final result for the characteristic speed is finite. 

We should emphasize that, by rescaling uq —>■ Iuq, the above calculation also provides the speed 5'outer 
in Eq. (I4.5t of section 0] More precisely, upon neglecting degenerate branches in Figs. I4.1h -c. the values 
of the Riemann invariants that determine ^^3^ in Fig. I3.2h are exactly those of the Riemann invariants that 
determine s'^ in case (i) (cf. Fig. 14. Ih ) and in cases (ii) and (iii) (cf. Figs. I4.1b .c') of section |4] upon 
replacing uq —>■ Iuq. 

Degenerate genus-2 calculations: Equations (I3.2> and (I4.8> . We now want to calculate in case (ii) 
of section|3j as shown in Fig. 13.2b . This is done by evaluating S5{ri,... ,r(,) = ^(rs), with s{r) still given by 
Eq. (lA.nt with g = 2 and with rg = — ri = uq + 2qo, = —r^, = —uq + 2<7o, = —r2 = r4 + £, and then 
taking again the limit £ — > 0+. In general, Eqs. (!A.14t yield the coefficients tti a2, Yi and 72 for the g = 2 
case as: 

a,-+i = (-l)^-+i^([3,7] - iai [2,7]) , (A.24a) 
Yy+i = (-1)''^' [T^ ([4.7] - [3,;] + \{o2- \oj^ [2,7]) , (A.24b) 



7 = 0, 1, where 



[/,;•] =det(^| ;•]. (A.24C) 



In our case, thanks to the symmetry of the initial datum for the invariants, Eqs. ( IA.13t give ai = and 
O2 = —rl — rl — rl = (— + 4mo^o — 12^q) + 2(mo — 2^o)£ — Upon making the substitutions 2 = ^3 + 
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(r4 — rs) sin^ 6 on and z = + (re — rs) sin^ 6 on /| and using the reflective symmetry of the initial datum, 
the integrals /["j, • • • 2 form: 

n/2 pns2<:9fi 

/f = (-l)M / , ^e, (A.25a) 



5 



r2cos2 2e)(r2-r2cos2 2e) 



, '^/2 (rs + Asin^e)*' 

/| = / ^ V a ; ^ ^ (A.25b) 

(rs - r4 + A sin^ 6) (r4 + rs + A sin^ 6) (2 rj + A sin^ 6) (rs + re + A sin^ 6) 





where A = rg — rs . Then, evaluating the above integrals, substituting Eq. (IA.25I) and Eq. (IA.24t into 
Eq. dA.nt and taking the hmit £ ^ 0+, one then obtains the desired result, namely Eq. ( I3.2I ). Again, 
note that even though some of the integrals are divergent in the limit £ ^ 0+, the result for the characteristic 
speeds is finite. 

Again, we emphasize that, b^ rescaling qo go + ^"o and uq —>■ Iuq, the above calculation also provides 
both of the speeds s^^l^^ and s[J^^^ in Eq. (I4.7t of section |4] More precisely, it is easy to see that, upon 



neglecting degenerate branches in Figs. I4.1b -d. the values of the Riemann invariants that determine in 
Fig. 13.2b are the same as those that determine in Fig. l4.1b -c qo —>■ qo + ^mq- Similarly, upon neglecting 
degenerate branches Figs. l4. Ib -d one sees that the values that determine s'^ in Fig. 13.2b also coincide with 
those that determine Sj in Fig. 14. lb and s'^ in Fig. 14. lb upon uq Iuq. 
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